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We construct exactly soluble lattice models for fractionalized, time reversal invariant electronic 
insulators in 2 and 3 dimensions. The low energy physics of these models is exactly equivalent to a 
non-interacting topological insulator built out of fractionally charged fermionic quasiparticles. We 
show that some of our models have protected edge modes (in 2D) and surface modes (in 3D), and 
are thus fractionalized analogues of topological insulators. We also find that some of the 2D models 
do not have protected edge modes - that is, the edge modes can be gapped out by appropriate time 
reversal invariant, charge conserving perturbations. (A similar state of affairs may also exist in 3D). 
We show that all of our models are topologically ordered, exhibiting fractional statistics as well as 
ground state degeneracy on a torus. In the 3D case, we find that the models exhibit a fractional 
magnetoelectric effect. 



I. INTRODUCTION 

One of the more surprising discoveries of the past 
decade has been that time-reversal invariant band insula- 
tors come in two kinds: topological insulators and trivial 
insulators. These two families of insulators exist in both 
twoi^— and three^"— dimensional systems. They are dis- 
tinguished by the fact that the interface of a topological 
insulator with the vacuum always carries a gapless edge 
mode (in two dimensions) or surface mode (in three di- 
mensions) , while no such protected boundary modes exist 
for the trivial insulator. 

Though much of our current understanding of topolog- 
ical insulators has focused on non-interacting or weakly 
interacting systems, it is natural to consider the fate 
of this physics in the presence of strong interactions. 
Strongly interacting insulators can be divided into two 
classes: systems that can be adiabatically connected 
to (non-interacting) band insulators without closing the 
bulk gap, and those that cannot. In the former case it has 
been shown (explicitly in 2D^, and implicitly in 3D^) that 
the gapless boundary modes of a topological insulator are 
stable to strong interactions as long as time reversal sym- 
metry and charge conservation are not broken (explicitly 
or spontaneously). Therefore there is a well-defined no- 
tion of interacting topological insulators in systems that 
are adiabatically connected to band insulators. 

Here we will consider the second possibility: strongly 
interacting, time reversal invariant electron systems 
whose ground state cannot be adiabatically connected 
to a band insulator. The same question can be posed: 
do some of these systems have protected gapless edge 
modes? This question is particularly interesting in light 
of the fact that such phases can be fractionalized, leading 
to a great diversity of possibilities. That is, such phases 
need not have excitations that resemble electrons; in gen- 
eral, the quasiparticles may have fractional charge and 
statistics. 

Our understanding of these fractionalized insulators 



is, however, limited. Focusing on the two dimensional 
case, Ref. [l^ analyzed a class of strongly interact- 
ing toy models^ where spin-up and spin-down electrons 
each form fractional quantum Hall states with oppo- 
site chiralitiesJ^ Ref. [13 concluded that some of these 
strongly interacting, time reversal invariant insulators 
have protected edge modes, while some do not. The 
two kinds of insulators were dubbed "fractional topolog- 
ical insulators" and "fractional trivial insulators" since 
they are analogous to non-interacting topological and 
trivial insulators, but they contain quasiparticle exci- 
tations with fractional charge and fractional statistics. 
These models demonstrate that fractionalized analogues 
of topological insulators are possible in principle, but 
they do not exhaust all the possibilities for these phases. 

In the three dimensional case, even less is known. Refs. 
[T2I and [1^ used a parton construction to build time rever- 
sal invariant insulators with a fractional magnetoelectric 
effect. However, this work did not prove that these states 
have protected surface modes (though Ref. did con- 
jecture that this is the case). Also, Refs. Il2l andlisl did 
not construct a microscopic Hamiltonian realizing these 
phases~a standard limitation of parton or slave particle 
approaches. 

In this paper, we address both of these issues. First, we 
construct a set of exactly soluble lattice electron models 
- in both two and three dimensions - that realize time 
reversal invariant insulators with fractionally charged ex- 
citations. Second, we prove that some of these fractional- 
ized electronic insulators have protected edge or surface 
modes (that is, we show that perturbations cannot gap 
out the boundary modes without breaking time reversal 
symmetry or charge conservation symmetry, explicitly or 
spontaneously). In this sense, these models provide con- 
crete examples of "fractional topological insulators" in 
both two and three dimensions. An added bonus of our 
analysis is that the argument we use to establish the ro- 
bustness of the edge or surface modes is not specific to 
our exactly soluble models, and can be applied equally 
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well to more general fractionalized (or unfractionalized) 
insulators. 

The low energy physics of our models is exactly equiv- 
alent to a non-interacting topological insulator built out 
of fractionally charged fermions. Surprisingly, however, 
some of the models do not have protected boundary 
modes. Specifically, we find that in some 2D models 
(namely those for which our protected-edge argument 
breaks down) the edge modes can be gapped out by 
adding an appropriate time reversal invariant, charge 
conserving perturbation. In the 3D case, our understand- 
ing is more limited: we can prove that some of the models 
have protected surface modes, but we do not know the 
fate of the surface modes in the remaining models. 

The models that we construct and solve are rather far 
from describing systems that are presently accessible to 
experiments. Nevertheless, they are of value for two rea- 
sons. First, they allow for a study of matters of principle, 
such as the existence of fractional topological insulators, 
their properties and their stability. And second, since 
the phases we consider are robust to arbitrary deforma- 
tions of the Hamiltonian that do not close the bulk gap 
(or break time reversal or charge conservation symme- 
try), these phases may also be realized by models that 
are significantly different from the ones discussed here. 

The paper is organized as follows. In section |lT] we 
describe the basic physical picture and summarize our 
results. In section IIIIl we construct models in the 2D 
case. In section HVl we analyze the physical properties of 
these models, including the structure of the edge modes 
and the topological order in the bulk. In section |Vl we 
consider the same models in the 3D case. We then an- 
alyze the physical properties of the 3D models in sec- 
tion I VII including the structure of the surface modes, 
the topological order in the bulk, and the nature of the 
magnetoelectric efi^ect. In the final part of the paper, 
section fVin we investigate whether the boundary modes 
in our models are robust to arbitrary time reversal in- 
variant, charge conserving perturbations. The Appendix 
contains some of the more technical calculations. Table 
H] lists the various symbols that we use in the text. 



Symbol 


Description 


Slprf inn 


Nsitc 
plaq 


number of sites s in the lattice 
number of links (ss') in the lattice 
number of plaquettes P in the lattice 


M 
rm 
rm 


m 

bl 
b^ , 

ss' 

ris 

Qs 
Bp 

as 

Uss' 

Qs 
bp 

\qs,bp) 

Qch 


integer parameter in boson model 
boson creation operator on site s 
boson creation operator on link (ss') 
boson occupation number on site s 
boson occupation number on link (ss') 
cluster charge on site s in boson model 
ring exchange term on plaquette P 
sublattice weighting factor for site s 
hopping term on link (ss') 
eigenvalue of Qs 
eigenvalue of Bp 

simultaneous eigenstate of Qs, Bp 
simultaneous eigenstate \qs,bp — 1) 
electric charge of charge excitation 




mm 


ITTT Al 
liii Al 






iTTTAl 


ITTT Al 

111! Al 




IlIIAl 


nil A 11 


IIIlAll 


IJII A 11 


nil A 21 


nil A 21 


ct 

Qs 

k 

qs 

\qs,bp,nsa} 
|?iso-, elec) 

In-sir) 

qt 

Hhop 


electron creation operator 
electron occupation number 
total number of electrons on site s 
cluster charge in electron model 
integer parameter in electron model 
eigenvalue of Qs 

simultaneous eigenstate of Qs, Bp, risa 
electron state with occupation {usa} 
eigenstate Ids —0,bp = l,nscr) 
electric charge of fermion excitation 
hopping term for fermion excitations 
hopping amplitudes for fermions 
creation operator for fermions 


ITTT F!l 
liii i3l 


InTBl 


ITTT F!l 
liii iJl 


ITTT P!l 
liii i3l 


InTBl 


HUB 11 


ITTT F! 1 1 

liii D il 


HUB 21 


HUB 21 


IIIIB2I 


ITiiq 


Imcl 


IlIICl 


dchfi 

6if,fl 


Mutual statistics of charges and fluxes 
Statistics of fermions and fluxes 


ITVBI 


llVBl 


e* 


smallest charged excitation 


IVICI 



TABLE I. List of symbols, their description, and the section 
where they are deflned. 



II. SUMMARY OF RESULTS 

This section is aimed at introducing the reader to our 
exactly soluble models, and to the main results we find by 
analyzing these models. We will emphasize the physical 
picture, leaving the detailed calculations to the following 
sections. 



A. Constructing exactly soluble models for 
fractional topological insulators 

To obtain candidate fractional topological insulators, 
we build models with two important properties: (1) frac- 
tionally charged fermionic quasiparticles and (2) a topo- 



logical insulator band structure for these excitations. 
Our construction has three steps. In the first step we 
construct lattice boson models with fractionally charged 
bosonic excitations. In the second step we add elec- 
trons to the lattice, and define an electron-boson inter- 
action that binds each electron to fractionally charged 
excitations of the bosonic model, thus creating a frac- 
tionally charged fermion. In the third step we construct 
a hopping term on the lattice that allows the fraction- 
ally charged fermion to hop between lattice sites without 
exciting other degrees of freedom. We then choose the 
hopping terms so that these fermions have a topological 
insulator band structure. 

The boson models we construct are similar in spirit 
to the "toric code" model and its Zm generalizationsii. 
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but are built out of bosonic charged particles whose to- 
tal charge is conserved. In this sense, these models are 
a hybrid between the toric code model (which is exactly 
soluble but not charge conserving) and the fractionalized 
bosonic insulators of Refs. [l5l - ll9l (which are charge con- 
serving, but not exactly soluble). 

The construction of the models is based on the follow- 
ing recipe. We consider a system of bosons that live on 
the sites s and links {ss') of a bipartite square (or in 
3D, cubic) lattice. We construct a bosonic Hamiltonian 
composed of two parts: Hb ^ Hi + H2- Each term has 
an associated energy scale whose magnitude is of minor 
significance to our discussion. Both, however, depend on 
an integer parameter m which plays a crucial role, as it 
determines the fractional charge carried by the quasipar- 
ticle excitations. 

The first term Hi is the "charging" Hamiltonian. This 
term depends only on the number of bosons on each site 
Us and on each link n^s'. Each boson is made of two 
electrons, of a charge e each. The charging Hamilto- 
nian assigns different energies to different charge config- 
urations {ns,rtss/}, by coupling the charge on a site to 
the charges on the four links neighboring the site. The 
spectrum of Hi is discrete, as expected from a charging 
Hamiltonian. The spectrum is also highly degenerate, 
since many charge configurations have the same energy 
cost. In fact, the number of degenerate eigenstates of 
the lowest eigenvalue of Hi is 771^""''"^='"'+^ with Mink 
being the number of links in the lattice and A^sitc being 
the number of sites. 

The second term H2 is the "hopping" Hamiltonian. 
This term makes bosons hop between neighboring lattice 
sites and links. A crucial aspect of our model is that the 
two parts are mutually commuting: [Hi^H^] = 0. Thus, 
the hopping Hamiltonian H2 only has matrix elements 
between degenerate states of the charging Hamiltonian 
Hi. and splits the degeneracy for the ground state. 

As we want to build an insulator, we need the ground 
state of Hb to be separated from the excited states by 
a finite energy gap. Furthermore, because we want frac- 
tionally charged excitations, Hb must be topologically 
ordered^Si^i (in gapped systems, fractional charge implies 
the existence of topological order) . The presence of topo- 
logical order means that the degeneracy of the ground 
state must depend on the topology of the systemi^"— 
More specifically, we need the degeneracy of the ground 
state to be independent of the system size, and to be 
different for a system with open and periodic boundary 
conditions. 

The first condition ~ existence of an energy gap - is 
guaranteed in our model by having the spectra of the 
charging Hamiltonian Hi and the hopping Hamiltonian 
H2 discrete. Note that this is not a common feature 
to hopping Hamiltonians. The continuous spectrum of 
the Josephson Hamiltonian is a representative example 
to the contrary. To make the spectrum discrete, we need 
to choose a carefully tailored hopping operator. While a 
conventional hopping Hamiltonian allows a single particle 



to hop between two neighboring sites, the hopping term 
we introduce allows only for a simultaneous correlated 
hopping of several particles around a single plaquette. 

The second condition - a ground state degeneracy that 
depends on the topology of the system - is a consequence 
of the way that the hopping Hamiltonian splits the de- 
generacy of the ground state of the charging Hamiltonian. 
For example, consider the case of the 2D system defined 
on a torus. Each of the terms in H2 describes hopping 
around one of the A^piaq plaquettes of the lattice and only 
one out of 77i^pi.->q~i ground states of Hi is also a ground 
state of H2- Thus, the ground state degeneracy of the 
Hamiltonian Hb on a torus is TTi^iink-A^sito-Wpiaq-i-s^ gy 

Euler's theorem, this number is exactly m^. A similar 
calculation in a 2D open geometry yields a ground state 
degeneracy of 1. In the 3D case, the analysis is similar. 
One finds that the ground state degeneracy in a 3D open 
geometry is again 1, while on a 3D torus it is nv^ . 

This counting agrees with the generalized "toric code" 
model with gauge group G ~ Z,n^ The quasiparticle 
excitations of the boson model are also similar to the 
toric code: there are two types of quasiparticle excita- 
tions - "charge" particles and "fiux" particles - which are 
individually bosons but have fractional mutual statistics. 
Also, like the toric code model, the boson model docs not 
have gapless edge modes. The main difference from the 
Zm toric code model is that the "charge" quasiparticles 
carry a fractional electric charge, 2e/m. 

After constructing the bosonic models, we next intro- 
duce single electron degrees of freedom that live on the 
lattice sites. The electrons couple to the bosons through 
the charging energy, and the electron-boson coupling is 
characterized by a second integer parameter k. We design 
this coupling so that it energetically binds an electron to 
a composite of k fractionally charged bosonic quasiparti- 
cles, each carrying charge 2e/m. The resulting composite 
particle then has a fractional charge oi qi ~ e{l + 2k /m), 
and follows fermionic statistics. We denote the Hamilto- 
nian of this modified lattice model by H^. 

In order for the composite particle to be a stable de- 
gree of freedom, it must be able to hop between lattice 
sites "in one piece" , i.e. without affecting the other types 
of excitations. In the final step of the construction, we 
find a hopping term ffhop that does just that. We then 
add Hhop to the Hamiltonian i/g, choosing the hopping 
amplitudes so that the composite particles have a band 
structure of a topological insulator. The energy gap be- 
tween the bands is a parameter of i?hop, and we assume 
it to be much smaller than the energy gap of the bosonic 
excitations. 

This construction results in a system of non-interacting 
fermions of spin-1/2 and charge qt = e(l -|- 2k /m) in 
a topological insulator band structure in either two or 
three dimensions. The smallest charged excitation in the 
system carries a charge e* = 2e/m when m is even, and 
a charge of e* = e/m when m is odd. In the former 
case, this is a bosonic excitation. In the latter, it is a 
composite of fermionic and bosonic excitations. 
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B. Properties of the models 

1. The two-dimensional case 

In two dimensions our models realize quantized spin 
Hall states, with a pair of gapless edge modes and a spin- 
Hall conductivity of + The topological or- 
der characterizing the states originates from the bosonic 
models underlying them. The ground state degeneracy 
on a torus is m?. In the bulk there are three types of 
excitations: the bosonic charge excitation with electric 
charge 2e/m, the bosonic flux excitation which is neutral, 
and the fcrmion excitation with charge qf = e(l + 2k/m). 
We find the flux excitation to have a non-trivial mutual 
statistics with the other two types of excitations. When 
a bosonic charge excitation of charge 2e/m winds around 
a flux excitation, it accumulates a phase of 27r/m. Con- 
sequently, when a fermion, which is a composite of an 
electron and k bosonic charge excitations, winds around 
a flux particle, it accumulates a phase of 27rfc/m. 

In certain limits the only active degrees of freedom are 
those of the fermions at the edge, where a gapless mode 
exists. The system can then be described as a topo- 
logical insulator built out of non-interacting fermions of 
fractional charge qt = e{l + 2k/m). In particular, this de- 
scription holds when the system is driven at low frequen- 
cies and long wave lengths by a weak electromagnetic 
field or when thcrmodynamical properties arc probed at 
low temperatures. Under these conditions, the system 
would show a two-terminal conductance of 2q^/h, the 
shot noise associated with tunneling between edges would 
correspond to a charge of gf , and the heat capacity would 
be linear in temperature and proportional to the system's 
circumference, as expected from a 2D topological insula- 
tor of non-interacting charge qf fermions. 

When deviating from these conditions, the bosonic de- 
grees of freedom can become active. Examples include 
the application of a magnetic fiux of the order of a flux 
quantum, $o — hc/e, per plaqucttc, the application of 
bias charges of order of e/2m to particular sites and the 
application of an electromagnetic field at frequencies that 
correspond to the gap to bosonic excitations. 

2. The three dimensional case 

In three dimensions our models are strong topological 
insulators built out of charge = e(l + 2k/m) fermions, 
with a gapless Dirac cone on each surface. When time 
reversal symmetry is broken on the surface, the models 
exhibit a surface Hall effect with a fractional Hall con- 
ductivity of q^/2h. As in the two dimensional case, the 
topological order in the 3D model originates from the 
topological order of the underlying bosonic system. The 
charged excitations carry electric charges of 2e/m and 
e{\-\-2k/m) and are identical to those in the 2D case, but 
the fiux excitation becomes a flux loop rather than the 
point particle it is in 2D. The ground state degeneracy 



on a 3D torus is . Again, in certain limits the bosonic 
degrees of freedom may be neglected and the only active 
degrees of freedom are the fermionic ones. The condi- 
tions for these limits to hold are similar to those of the 
two dimensional case. 

The bosonic degrees of freedom are active in several 
cases, one of which is of particular interest. In a 3D topo- 
logical insulator of non-interacting electrons, a magnetic 
monopole in the bulk of the insulator binds a half integer 
electric charge.—"— Hence, a monopolc/anti- monopole 
pair - which may be created by a finite-length solenoid 
carrying a fiux quantum $o a-nd positioned within the 
bulk - creates an electric dipole with a half-integer elec- 
tric charge at its ends. In our model we find that such 
a solenoid leads to a dipole with a charge which is a 
half-integer multiple of q^/e. Unlike the non-interacting 
case, however, the energy involved in creating the dipole 
is proportional to its length - indicating that the two 
ends of the dipole cannot be effectively separated from 
one another. The two ends of the dipole can be sepa- 
rated only when the fiux carried by the solenoid is e/e* 
fiux quanta. Furthermore, because we could presumably 
trap any number of additional charge e* quasiparticles 
near the ends of the solenoid by adding an appropriate 
local potential, the only quantity which is independent 
of microscopic details is the monopole charge modulo e*. 
Calculating this quantity, we find that the charge at the 
end of the (e/e*)3>o solenoid is a half- integer multiple 
of e* for the models where qt/e* is odd, and an integer 
multiple of e* for the models where qt/e* is even. 



C. The stability of the edge or surface modes 

In conventional topological insulators, the edge or sur- 
face modes are protected as long as time reversal symme- 
try and charge conservation are not brokeni^ If either of 
these symmetries is broken, e.g. by a Zeeman magnetic 
field that couples to the electron spin or by a proximity- 
coupling to a superconductor that allows for Cooper pairs 
to tunnel into and out of the edge or surface modes, these 
modes may be gapped. The breaking of time reversal 
symmetry may be spontaneous rather than explicit, in- 
duced for example by the Fock term of electron-electron 
interaction. The stability of the edge or surface modes 
to perturbations that do not break these symmetries is 
the distinguishing feature of topological insulators in 2D 
and strong topological insulators in 3D. 

An important question is whether the phases we study 
here have protected edge or surface modes similar to con- 
ventional topological insulators. We find that some of 
the models do indeed have edge or surface modes pro- 
tected by time reversal symmetry and charge conserva- 
tion, while some do not. (Independent of this difference, 
all the models arc topologically ordered, as demonstrated 
by their topological ground state degeneracy). 
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1. The two dimensional case 

In the 2D case, wc find that our models conform to the 
general rule derived in Ref. that is, the edge modes 
are protected if and only if the ratio (JsH / e* is odd, where 
UsH is the spin-Hall conductivity in units of e/27r and e* 
is the elementary charge in units of e. In our models, 
this criterion is equivalent to the condition that the ratio 
Qf/e* is odd. 

We establish the stability of the edge modes for the 
models with odd qt/e* by a general flux insertion argu- 
ment similar to the used in Ref. ITol and establish the in- 
stability in the case of even qt/e* by explicitly construct- 
ing the perturbations whose combination gaps the edge. 
This combination is rather interesting. As defined, the 
models have two fermionic edge modes of opposite chiral- 
ities - the bosonic excitations are gapped at the edge. In 
order to gap the fermionic edge modes, we introduce one 
perturbation whose role is to close the gap of the bosonic 
excitations at the edge, and then two additional pertur- 
bations that couple the bosonic and fermionic modes, 
gapping them both. 

For the closure of the bosonic gap at the edge we apply 
a perturbation aimed at turning the edge of the bosonic 
system from an insulator into a superfluid. The natural 
way of doing that is by introducing a hopping Hamilto- 
nian that allows fractionally charged bosonic excitations 
at the edge to hop from one site to another. When the 
hopping term is strong enough it can overcome the charg- 
ing term described by Hi , thereby closing the gap at the 
edge. As for the perturbations that couple the bosons 
and the fermions at the edge, the first such perturbation 
breaks a spinless boson of charge 2e into two electrons 
of opposite spin directions on the same lattice site. The 
second of these perturbations flips the direction of an 
integer number of electrons' spins, while simultaneously 
operating on the flux degrees of freedom on the edge. 
Both of these perturbations make use of the bosonic de- 
grees of freedom and therefore do not have analogues in 
non-interacting electron systems. 



2. The three dimensional case 

Just as in the 2D case, we find that the 3D models with 
odd Qf/e* have protected surface modes. We establish 
this result using a 3D generalization of the flux insertion 
argument of Ref. [l^. We note that this argument is 
of interest beyond the particular models discussed here, 
and can be applied to more general fractionalized and 
conventional insulators. Unlike the 2D case, we are not 
able to determine the stability of the surface modes for 
the models with even q{/e*. Addressing this question 
requires either the construction of specific perturbations 
that gap out the surface, or an argument proving that 
the surface modes are protected. 




s <ss'> s' 



FIG. 1. In the lattice boson model, bosons live on both the 
sites s and links {ss') of the square lattice. The Hamiltonian 
Hb m is a sum of a Qs term ((2]), which acts on four links 
{ss') and one site s, and a Bp term which acts on the 
four sites and links adjacent to a plaquette P. The Bp term 
is a product of four link operators Uggi (O which each act on 
the sites s,s' and the link {ss'). 



III. LATTICE MODELS FOR 2D FRACTIONAL 
TOPOLOGICAL INSULATORS 

A. Step 1: 2D lattice boson models with fractional 
charge 

In this section we describe a collection of exactly 
soluble lattice boson models with fractionally charged 
excitations — one for each integer m > 2. The models 
can be defined on any bipartite lattice in 2 or higher di- 
mensions. Here, for simplicity, we will focus on the case 
of the square lattice. Later, when we construct 3D mod- 
els, we will consider the cubic lattice case. 

The basic degrees of freedom in these models are charge 
2e spinless bosons which live on the sites s and links (ss') 
of the square lattice. We denote the boson creation oper- 
ators on the sites and links by 6j and b],^, and the corre- 
sponding boson occupation numbers by and n^s' . The 
Hamiltonian Hb can be written as a sum of two terms, 
one associated with sites s, and the other associated with 
plaquettes P of the square lattice (Fig. [T]): 

Hb = Hi + H2 

= vY.Q"^^Y.(Bp + bU (1) 

s P 

We will take u,V > but otherwise arbitrary. The Qs 
term is a "cluster charge" term which measures the total 
charge on the site s and the four neighboring links (ss') 
with appropriate weighting factors. It is defined as the 
sum 

Qs = as ^ Uss' + m - Us (2) 

s' 

where 

as^l' ^'^^^ (3) 

I m — 1 if s G i? 
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and A and B are the two sublattices of the square lat- 
tice. Since V is positive, VQl describes a short range 
repulsive interaction between the bosons. This interac- 
tion breaks the sublattice symmetry between the A and 
B sublattices, except in the case m = 2. The Bp term 
can be thought of as a ring exchange term. It is defined 
as the product 

Bp = U12U23U34U41 (4) 
where Uss' is a boson hopping term on the link {ss'): 



(5) 



The hopping term Uss' describes processes where bosons 
hop from the sites s, s' to the link {ss') and vice versa. 
It is designed so that it has two special properties. First, 
Uss' changes the number of bosons on the site at the 
center of the link {ss') by ±1 (mod m) with the -I- sign 
when s € B and the — sign when s € A. This change is 
compensated by a corresponding increase or decrease in 
the number of bosons in the two neighboring sites s, s' 
so that the total number of bosons is conserved. Second, 
Uss' decreases the cluster charge Qs by 1 and increases 
the cluster charge Qs' by 1 and doesn't affect the charge 
on any other site: 



[Qr, Uss'] = {Srs' - Srs)Us. 



(6) 



An important consequence of this relation is that Qs 
commutes with the product of Uss' around any set of 
closed loops, and in particular, 



[Qs,Bt 







(7) 



Equation ([7]) is at the root of why our system is an in- 
sulator: the Bp operator has no effect on the cluster 
charges Qs and hence does not provide for the long- 
distance transport of electric charge. 

The final component of the model is our definition of 
the boson creation operators bl , bl^, . For the site bosons 
&|, we use a rotor representation, letting &| = e*^° with 
[6s-, ns] = i. The boson occupation number on the sites 
can therefore be any integer, Us S (—00,00). On the 
other hand, we take the link bosons bl^, to be a kind of 
generalized hard-core boson, restricting the boson occu- 
pation number to Uss' G {0,1,..., to — 1}, and defining 
^Ic' to be the m X m matrix 



/O 1 
1 


VO 



0\ 


1 

0/ 



(8) 



when written in the normalized number basis {|to — 
1),...,|0)} on the link (ss'). We note that while these 
generalized hard-core bosons are unconventional, they 



can arise as an effective description of a conventional 
boson system in an appropriate limit. For example, if 
the number of bosons on the link {ss') is very large but 
is restricted to a set of m contiguous values {TV, Af + 
1, ...,J\f+m—l} by appropriate energetics, then the above 
definition of 6^^, becomes a good approximation to con- 
ventional bosons (up to an overall normalization factor). 
To summarize, the full Hilbcrt space for our model is 
spanned by the occupation number states \ns, nss') with 
Us G (—00, 00) and < Uss' < m — I. 



1. Solving the boson model 

We will now show that the Hamiltonian ([T|) is exactly 
soluble, and compute its exact energy spectrum. We are 
already part way there, having established that Qs,Bp 
commute with each other ([7]). Next, we note that 

[Uss',Urr']=0 , Ul^U;;^Us's (9) 

from which it follows that 

[Bp,Bp,]^[Bp,BU = . (10) 
Combining these results with the obvious relation 



[Qs,Qs 



0, 



(11) 



we conclude that {Qs,Bp,Bp} all commute, and there- 
fore can be diagonalized simultaneously. 

The simultaneous eigenstates of these operators can be 
labeled as \qs,bp), where 



Qs\qs,bp) = QsIqs, bp) 
Bp\qs,bp) = bp\qs,bp) 
Bl,\qs,bp) = b*p\qs,bp) 

The corresponding energies are 

E-VY^q^-^Y.(f^P + b*p) 



(12) 



(13) 



It is clear from the definition ^ that Qs has integer 
eigenvalues so qs is an integer. Similarly, using the 
fact that Uss' changes the occupation number Uss' by 
±1 (mod m), we can show that 



= 1 



(14) 



so bp must be a mth root of unity. This relation guaran- 
tees what wc promised in section |TT) the spectrum of H2 
is discrete. 

The only remaining question is to determine the degen- 
eracy of each of the qs,bp eigenspaces. This degeneracy 
depends on the geometry we consider. We first consider 
the case of a rectangular piece of square lattice with open 
boundary conditions (Fig. [HJa)). We show in Appendix 
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FIG. 2. We consider the lattice boson model Hb (O in two 
geometries: (a) a rectangular geometry with open boundary 
conditions, and (b) a periodic (torus) geometry. In the rect- 
angular geometry (a), the Qs operators at the corners act on 
two links (ss'), while those at the edge act on three links. 



I A II that in this geometry, there is a unique eigenstate for 
each collection of {qs, bp} satisfying the global constraint 



(mod m) 



(15) 



In other words, {qs,bp} are independent and complete 
quantum numbers, except for this single global con- 
straint. A similar result holds for a (periodic) torus geom- 
etry (Fig. IHb)). In this case, we find that there are rn? 
states for each collection of {qs, bp} satisfying ([T5|) as well 
as the additional constraint Jlp^f ~ 1 i^^'^ Appendix 



72]). This degeneracy is a consequence of the fact that 
the system is topologically ordered (see Section I IV Bp . 
Below we will focus exclusively on the open boundary 
condition geometry, unless otherwise indicated. 

Putting this all together, we conclude that the ground 
state of dD) is the unique state with Qs = 0, 6p = -1-1 ev- 
erywhere. There are two types of elementary excitations: 
"charge" excitations where qs = I for some site s, and 
"flux" excitations where bp = g^'^^f™ for some plaquette 
P. The total number of charge excitations, in the bulk 
and edge together, must sum up to modulo m. The 
charge excitations cost an energy of V, while the flux 
loop excitations cost an energy of u{l — cos(27r/m)). In 
particular, as long as u,V > 0, then the ground state is 
gapped. 



2. Fractional charge in the boson model 

An important property of the boson model ([1]) is that 
the charge excitations carry fractional electric charge 
Qch = 2e/m. We can derive this result by explicitly cal- 
culating the electric charge distribution in these states. 
Consider an eigenstate \qs) = \qs,bp = 1) with some ar- 
bitrary configuration of charges {qs}, and with no fiuxcs. 
It is straightforward to show that the (un-normalizcd) 
microscopic wave function for this state in the occupa- 
tion number basis is given by 



{nss',ns\qs) = 



otherwise 



-|- mrLs = qs for all s 



(16) 



(Again, to be precise, this wave function applies to a sys- 
tem with open boundary conditions) . This wave function 
has a nice property: if we were to measure the occupa- 
tion number riss' on any link (ss'), then we would find 
each of the m possible values, riss' = 0, 1, ...m — 1, with 
equal probability, 1/m. This is true on every link ss', in- 
dependent of the configuration of charges {qs}- It follows 
that the expectation value of Uss' in the state \qs) is 



{'^ss'}qs 



O + 1 + ... + m- l 



m — 1 



(17) 



Using this result, together with the constraint 
as Uss' + mus = qs, wc deduce 



m 

zas{m — l)/2 



(18) 



where z is the coordination number of lattice (z = 4 for 
the square lattice). 

Wc are now in a position to compute the fractional 
charge q^h. Consider the case of an isolated charge 
excitation — that is, a state where qs = 1 at some site 
So and = at all nearby sites. Denote this state by 
\qso = 1). Similarly, consider an eigenstate \qso = 0) 
where qs = both at sq and at all nearby sites. Then 
the electric charge qch carried by the excitation is given 
by the difference in expectation values 



qch = 2e[{Ns)i-{Ns)o] 



(19) 



where is an operator which measures the total num- 
ber of bosons in some large area S containing sq: 



Ns 



s.s' 



(20) 



Using the above results (|17m8p we derive 

{nss')i - {nss')o = 

{ns)i - {ns)o - ^ (21) 
m 

implying that g^h = Furthermore, we can see from 
these expressions that this charge is perfectly localized 
to the site sq. This perfect localization is specific to the 
exactly soluble model: in a generic gapped system we 
expect excitations to have a finite size of order the cor- 
relation length. 

Alternatively, we can derive the fractional charge using 
a simple identity: for any set of sites S in the square 
lattice, we have the relation 



Qs = m \'^ns+ ^ Us 
ses \ses s.s'es 

= mNs + ^ asUss' 

sG5,s'e5= 



E 



UsUs 



ses.s'eS" 



(22) 
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Taking expectation values of both sides in the two states 
lis = 1): \Qs = 0)i a-nd subtracting gives 

{Ns)i-{Ns)o = -- V a,[{nss')i~{nss')o] (23) 
m ^ — ' 

To complete the calculation, we note that the second 
term on the right hand side vanishes in the limit that 
S becomes large, since in that case, the sum only in- 
volves links ss' that are far from sq, and the excess charge 
{{nss')i — {nss')a) must vanish at large distances from sq. 
It follows that Qch = 2e/m, as claimed. 

B. Step 2: 2D lattice electron models with 
fractional charge 

We are now ready to construct a model with fraction- 
ally charged spin-1/2 fermionic quasiparticle excitations. 
We accomplish this task by modifying the lattice boson 
model JT]) and coupling the bosons to additional (un- 
paired) electron degrees of freedom. 

The Hilbert space for this modified model is very simi- 
lar to the boson model, except that we introduce an addi- 
tional spin-1/2 electron degree of freedom at each lattice 
site s. In addition, we now think of the lattice bosons 
as being charge 2e, spin-singlet pairs of electrons. This 
microscopic picture for the bosons is important conceptu- 
ally because ultimately we want a model for a fractional 
topological insulator which is constructed out of electron 
degrees of freedom. 

We will denote the creation and annihilation operator 
for the (unpaired) electrons by cl^,Csa, and their occu- 
pation number by Usa- We will use Us^e to denote the 
total number of the (unpaired) electrons on site s: 

Later, we will be interested in models with multiple or- 
bitals on each lattice site s. In that case, we will let a 
include both the spin and orbital degrees of freedom. 

In this notation, the Hamiltonian He for the electron 
model is a sum of three terms: 

He = vY,Ql-^J2(Bp + Bl.)-^iJ2ns. (25) 

s p sa 

where 

Qs ^ Qs - k ■ ns^e (26) 

and Bp,Qs are defined as in Eq. ([JIE]). We will take 
V,u> 0, but will consider both positive and negative /i, 
and arbitrary integer k. 

1. Solving the electron model 

The electron model ([25|) can be solved in the same way 
as the original boson model. Just as before, 

[Q„ Qs'] - [Bp, Bp,] = [Bp,Bl,] = [Q„ Bp] = (27) 



Also, it is clear that 

\Qs',ns„] = [Bp,nsa]^Q (28) 

Hence, we can simultaneously diagonalize 

{(3s,i?p,i?p,nso-}. Let \qs,bp,nsc) denote the si- 
multaneous eigenstates, where qs is an integer, 6p is a 
mth root of unity, and risa = 0, 1. The corresponding 
energies are: 

^ = ^ E 91 " i E(^^ + ^p) - E (29) 

s P sa 

By our analysis of the lattice boson model, we know that 
there is a unique state for each choice of qs,bp,ns^ sat- 
isfying the global constraint 

qs + k ^ User = (mod m) (30) 

S S(T 

(Again, we are assuming a geometry with open boundary 
conditions). 

Putting this all together, and taking /i < for sim- 
plicity, we conclude that the ground state is the unique 
state with ijs = 0,bp = l,nsa- = 0. The system has 
three types of elementary excitations: "charge" excita- 
tions with qs = 1 for some site s, "flux" excitations 
where bp = ^'"^1™ on some plaquette p, and spin-1/2 
"fermion" excitations, where = 1 for some site s and 
spin a =t, |. The charge excitations cost an energy of V , 
the flux excitations cost an energy of w(l — cos(27r/m)), 
and the fermion excitations cost energy — /i. In particu- 
lar, as long as w, y > and /i < 0, the ground state is 
gapped. 

2. Fractional charge m the electron model 

Our next task is to show that the fermion excitations 
carry fractional charge. To do this, we first need to intro- 
duce some notation. Let \ns„,elec) denote the electron 
occupation basis state 

|n..,dec)=n(cLnO) (31) 

s 

where |0) is the empty state. Also, let \nss' ,ns) denote 
the boson occupation state defined in section IIII Al A 
complete basis for the Hilbert space of the electron model 
is given by tensor product states \nsa,Glec) ® |ns,?^ss')- 
In addition to these general basis states, we will also 
find it useful to think about the set of eigenstates \nscr) = 
\qs = 0,bp = l,nsa-) made up of some arbitrary config- 
uration of fermions {usa} with no charge or flux exci- 
tations. (To construct these states we must impose the 
global constraint X^sa ""Sf^ = ^ (mod m)). The micro- 
scopic wave function for these states is given by 

\nscr) = \nscr,elec) iSi \qs = kns,e) (32) 
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where \qs) are the boson eigenstates defined in ^TE\\ . 

In order to compute the fractional charge carried by 
the fermion, it suffices to consider a state with an isolated 
fermion excitation — that is, suppose = 1 at some site 
sq and Usa- vanishes at all nearby sites. Denote this state 
by |"so(T = !)• According to the definition (|32l). this state 
is a tensor product of an electron state and a boson state: 

\nsoa- = 1) = \ns„a = 1, elec) (g) = k) (33) 

We can see that — 1, elec) consists of a single spin- 

a electron at site sq, while \qsg = k) corresponds to k 
"charge" excitations at site sq- Thus, the fermion exci- 
tation is a composite particle made of an electron and 
k charge excitations. To compute the total charge of 
the fermion, we need to add together the contributions 
coming from these two pieces. By our analysis of the 
fractional charge in the bosonic model, we know that 
each charge excitation carries charge 2e/m. On the other 
hand, the electron clearly has charge e. Adding together 
these two contributions, we conclude that the fermion 
excitation has charge 

9f = e(l + 2fc/m) (34) 



3. Time reversal symmetry and the electron model 

To construct candidate fractional topological insula- 
tors, it will be important to understand how the fcrmionic 
excitations in our model transform under time reversal. 
We use the usual convention for T, where the electron 
creation operators transform according to: 

r-- cl,^cl,cl^-el^ (35) 

In this convention, the bosons in (|25p transform trivially, 
since they are spin singlet pairs of electrons: 

bl, ^ bl, ,bl^bl (36) 

Applying these transformation laws, we see that the 
fermion excitations transform like spin-1/2 electrons. 



C. Step 3: Building candidate 2D fractional 
topological insulators 

In the last two sections, we have shown that the 
fcrmionic excitations of the electron model (j25p carry 
spin-1/2 and transform under time reversal just like elec- 
trons. In fact, they arc virtually indistinguishable from 
electrons except for the fact that they carry fractional 
charge qt ([M)) . Given these properties, it is easy to build a 
candidate fractional topological insulator: we simply put 
the fractionally charged fermions into a non-interacting 
topological insulator band structure. We accomplish this 
by adding a new term to the electron Hamiltonian 



(EH: 

H = He + ffhop 

\ s P sa ) 

+ i^hop (37) 

where 

H"hop = - ^ (t.ss'aa'c\,„,C,,\J^^^, -f h.c) (38) 

and Qs, Bp,Uss' are defined as before. The new term 
i/iiop gives an amplitude for the fermion excitations to 
hop from site to site without affecting any of the other 
degrees of freedom, as wc now show. Wc will assume 
that tss'tjtj' <^ u,V so that the bandwidth of the fermion 
excitations is much smaller than the gap to the bosonic 
excitations. 

We can understand the effect of -ffhop by computing 
the matrix elements of this operator between different 
eigenstates of , 

{q's^bp,n'^jHi,op\qs,bp,nsa) (39) 

This computation is considerably simplified by the fact 
that 

[Q.,ifhop] = [Sp,ifhop] =0 (40) 

implying that the matrix elements are only nonzero when 
q'g = (js, and b'p = bp. In what follows, wc specialize to 
the (js = q's = 0,bp = b'p = 1 case, since these are 
the lowest energy states and this is all we will need to 
understand the low energy physics. These states contain 
only fermions and no other excitations. As in (|32p . we 
will denote a state with some arbitrary configuration of 
fermions {rirr} using the abbreviated notation In^r) ^ 
\qs = 0, 5p = 1, Tirr) (Here r labels the sites of the lattice, 
while T =ti i labels the two possible spin states). To find 
the matrix elements (n^^lffhopl'TTr), we write 

cl,„CsaU^^,\nrT) = cl,^,Csa\nrT,elec) (E) U!:^,\qr = krir^e) 

(41) 

and then analyze each of these two pieces in turn. Using 
the explicit form of \qr) (fTB|) . we find 

Ut\qr) = \q:.) (42) 

where 

9r = Qr - kSrs + kSrs' (43) 

Similarly, wc have 

c\,„,Caa\nrT,dec) = ±\n'^^,elec) (44) 

where 

n'j.^ = Hrr - SrsSar + Srs'Sa'T (45) 
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and the ± sign depends on the ordering of the electron 
creation and annihilation operators in Eq. (j3ip . Com- 
bining these two results with (|4T|) . we derive 

cl,„Cs^U'^,,\nrr) = ±K^,e;ec) ® Ig^) = ±\n'^^) (46) 



where nj,^ is defined as in psj) . This relation establishes 
what we promised earlier: iJhop gives an amplitude for 
the fermion excitations to hop from site to site, but does 
not affect the other types of excitations. (We should not 
take these properties for granted: for example, if we had 
not included the operator U^^, in the definition of iJhop 
([38)) then iJhop would have affected the bosonic charge 
excitations) . 

Denoting the creation operators for the fermion excita- 
tions by , we conclude that the matrix elements of H 
within the low energy qs = 0,bp ~ 1 subspace are given 
by the free fermion Hamiltonian 

-ffoff = - ^ {tss'crcr'dl,^,dsa + h.C.) ~ fl'^dl^dscr (47) 

{ss') sa 

To complete the construction, we choose the hopping am- 
plitudes tss'aa' and the chemical potential ^ so that _ffoff 
is a non-interacting topological insulator. The low en- 
ergy physics is then described by a topological insulator 
built out of fractionally charged fermions. 

There are many possible choices for tg^'aa' , tJ-, but to 
be concrete we will focus our discussion on the following 
tight binding model on the square latticoS^. We con- 
sider a model with two orbitals, whose hopping matrix 
elements arc related by time reversal symmetry: letting 
(cr, cr') denote the fermion spin, and (1,2) denote the or- 
bital index, we have: 



,* 

''ss'aa' ,2 



(48) 



and there is no hopping matrix element connecting the 
two orbitals. For each spin, the hopping matrix elements 
for orbital 1 are: 

tss'crcr',! = ^ ^s-s',ei {td z + i\(Ji) + l5s-s',oi(K " z 
ei—x^y 

H = (49) 

where cr^^y^z are Pauli matrices acting on the spin in- 
dices. Here t parametrizes the spin-independent hopping 
terms, and A is a spin-orbit coupling. The constant k 
sets the band gap at the momenta (0, 0), (±7r, 0), (0, ±7r); 
in the regime < k < 2, the model is in a topologically 
insulating phase. 

In addition to being a topological insulator, the model 
P5|l conserves the total z component of the spin. The 
ground state consists of a filled spin-up band with Chern 
number v = +1, and a filled spin-down band with Chern 
number = — 1. As a result, the model exhibits a spin- 
Hall conductivity of 



(50) 



D. Effect of an electromagnetic field 

In this section, we investigate the response of the frac- 
tionalized insulator ([37]) to an applied electromagnetic 
field. We show that for weak, slowly varying fields, the 
model behaves like a non-interacting system of fraction- 
ally charged fermions. On the other hand, for stronger 
fields, we find that other, non-fermionic, degrees of free- 
dom contribute to the response. 

The first step is to understand how to incorporate a 
vector potential A into the Hamiltonian. As the model 
([37| contains charged bosons that live on links of the 
lattice in addition to those that live on the sites, we de- 
fine a lattice vector potential Ass',i, Ass', 2 for each of the 
two halves of each link {ss') (see Fig. [3]). The hopping 
operator Ugs' is then^l 



Uss' = (5t)"=-'6l,6^;,e2«e(l-a,)Aig2.eA. 



(51) 



~by-% bl, 



g2ie(l-Q,,)A2g2ieAi 



We can simplify this expression with the help of the uni- 
tary transformation 



Wa — exp 



2ie 
m 



Y^ris 

{ss') 



{asAss',1 - as'Ass',2) 



A little algebra shows that 



WaUss'W: 



bl) 



+ by-\ 



bl,b° 



J — A, 



(52) 



where Ass' = ^ss'.i + ^ss'.2 is the total vector potential 
on the link {ss'). In retrospect, this expression is to be 
expected as Uss' hops a charge 2e/m from s to s', and as 
such, should be multiplied by a phase factor e^**^^""' 
in the presence of an electromagnetic vector potential. 
Substituting this expression into H ([57]) . we find that 
the Hamiltonian can be written as 



s P 



where (j)p 
Hhop = 



A12 + A2 



A 



34 



hop 



AaI: 



h.C.) 

(53) 



its..'..'e''''-^-'cl,^,Cs.U^s' + h.C.) (54) 

{ss') 



and Qs, Bp,Uss' are defined as in the A = case. 

We now analyze this Hamiltonian in several cases. 
First, we consider the case where A is time independent 
and the fiux (jjp through each plaquette is small com- 
pared with a unit flux quantum. We proceed by simul- 
taneously diagonalizing {Qs, Bp, Bp,nsa}, denoting the 
eigenstates by \qs, bp, Usa). We then define rotated states 

\qs,bp,nsa,Ass') = WA\qs,bp,nsa) (55) 
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FIG. 3. In order to include an electromagnetic vector poten- 
tial in H (|37|) . we need to define a vector potential for each 
of the two halves of each link (ss'). We denote these vector 
potentials by A^s' ,i, A^^' ,2- 



We note that in the absence of the fermion hopping term 
iJhop, these states are exact eigenstates of H with ener- 
gies 



(56) 



When the flux (f)p is smah (specifically, \(t>p\ < $o/4), 
the lowest energy states are those with Qs = Q,bp = 1. 
We will denote these (purely fermionic) states using the 
abbreviated notation 



,Ass') = WA\qs = 0,6p = l^Urr) 



(57) 



To obtain the low energy effective Hamiltonian ffoff, we 
project H (|53)) to the subspace spanned by these states 
(|57p . By the calculation in the previous section, the ma- 
trix elements of (|53p within this subspace are given by 

H,B = - ^(t..'..'e^«'-^"'dL'rf- + h.c.) 

{ss'} 

-M^^dLrf- (58) 

This result proves that the model behaves like a non- 
interacting system of fractionally charged fermions, even 
in the presence of a weak, time independent vector po- 
tential A. 

In fact, the low energy effective Hamiltonian (|58p is 
also valid for a weak, time dependent A, as long as A 
varies slowly compared with the energy gap u,V of the 
bosonic excitations. One way to derive this result is to 
note that when A is slowly varying, we can make an adi- 
abatic approximation and can assume that the system 
always remains in the instantaneous low energy subspace 
spanned by the fermion states jrirr, ^ss')- The time evo- 
lution is therefore described by projecting the Hamilto- 
nian (|53p to the instantaneous low energy subspace, and 
wc again obtain the low energy effective Hamiltonian _ffcff 
([58| . As for Berry phase effects, one can check that 
non-abelian Berry connection ^ss' I^^^I'^'-'t', ^ss') 
reduces to a overall c-number phase factor and hence can 
be neglected for our purposes. (In other words. Berry 
phase effects only contribute a global phase factor to the 
time evolution of the wave function). 



We emphasize, however, that the above low energy ef- 
fective theory ([55)) is only valid when the flux through 
each plaquette is small compared to $o- Examining ([56l) . 
we see that when the flux becomes comparable to $0j 
the states with bp ^ 1 become energetically favorable. 
Hence, in this regime, we cannot describe the low energy 
physics in terms of the fermions alone: we also need to 
keep track of the other types of excitations. 

One case where these additional degrees of freedom 
are particularly relevant - and which will play an im- 
portant role in our analysis in section IVIII - is in flux 
insertion thought experiments. Imagine we take a ge- 
ometry with open boundary conditions, and we adia- 
batically increase the flux through a plaquette Pq from 
(/)Pq = to = $0i while keeping the flux through 
all other plaquettes constant at (j)p ~ 0. For simplicity, 
let us assume that the system is initialized in the state 
\qs = 0, 6p = l,f^so- = 0,Ass') when (ppg = 0, and let 
us neglect the fermion hopping term i/hop- From (|56p . 
we can see that there are two level crossings during the 
flux insertion process - one at (/)p„ = $o/4 and one at 
(f)Pg ~ 3$o/4. At the flrst level crossing, the state with 
bp„ = e-2iri/?ri becomes lower in energy than the state 
with bp„ = 1, while at the second level crossing, the 
state with bp„ = g^^'^*/™ becomes lower in energy than 
bpg = e^^TTj/m^ Pqj, ^ f^niiQ system size, we expect that 
these crossings will be avoided crossings, at least if we add 
a generic local perturbation to the Hamiltonian. There- 
fore, if we insert flux sufficiently slowly, the system will 
follow the lowest energy state at all times, finally evolv- 
ing into a state with bpg = e~^'^*/™ and (jypg = $o- One 
can check using the definition ([55]) that this final state 
is the same as the initial state. Hence, the insertion of 
a unit fiux quantum returns us to our original starting 
point - as expected from general considerations (e.g. the 
Byers-Yang theorem). 

An important point, however, is that the gap at these 
avoided level crossings vanishes in the limit that the pla- 
quette Pq is far from the edge of the system, since the 
states involved in the level crossings have different val- 
ues of 6pq and therefore only couple to one another via a 
process where a fiux quasiparticle tunnels from the edge 
to Pq. As a result, the above picture is only valid at ex- 
tremely long time scales. If we insert the fiux at a rate 
which is fast compared with the gap at the level crossing 
(but still slow compared with the bulk gap) the outcome 
is different. In this case, the system passes through the 
level crossings unaffected, leading to a final state with 
6po = 1 at $Pj, = $0- One can check using ([55]) that this 
final state corresponds to having 2 fiux quasiparticlcs on 
the plaquette Pq. This result shows that other degrees of 
freedom beyond the fermions come into play in fiux in- 
sertion experiments. In addition, it implies that we need 
to insert m fiux quanta (if m is odd) or m/2 fiux quanta 
(if m is even) for the system to return to its original con- 
figuration. In this sense, our models have a reduced flux 
periodicity - like fractional quantum Hall statesi^ 
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IV. PHYSICAL PROPERTIES OF THE 2D 
MODELS 

In the previous section, we introduced a Hamiltonian 
H ((571) whose low energy physics is exactly described by 
a non-interacting 2D topological insulator of fractionally 
charged fermions. We now describe the physical proper- 
ties of this fractionalizcd insulator, and derive an effective 
Chern-Simons field theory which summarizes them. 



A. Edge states 

One of the distinguishing features of non-interacting 
2D topological insulators is that they have gaplcss edge 
states which are protected by time reversal symmetry and 
charge conservatioujii^ii Here we discuss the analogous 
edge states for the 2D fractionalizcd insulator H ([57)) . 

Fortunately, the exact mapping between the low en- 
ergy physics of H and the non-interacting topological 
insulator H^s (HZ]) holds for any lattice, including ge- 
ometries with an edge. Hence, our problem reduces to 
understanding the edge structure of the insulator i?cff- 
We focus on the specific band structure (|49p , for simplic- 
ity. This insulator has two filled bands: a spin-up band 
with Chern number i/ = 1, and a spin-down band with 
Chern number v — Therefore, by the usual bulk- 
boundary correspondence, -ffcff must have one spin-up 
edge mode and one spin-down edge mode with opposite 
chiralitiesj^ Translating this result over to the fractional- 
ized insulator H , we conclude that H also has a pair of 
counterpropagating free fermion edge modes. The only 
difference from the non-interacting case is that the low 
energy fermions carry charge instead of charge e. The 
low energy edge Hamiltonian for H is thus of the form 

-ffcdgc = vd\{idx + qiA^)d^ - vd\{idx + qiA^x)di (59) 

where v is the velocity of the two edge modes. We leave 
the discussion of the stability of these edge modes to 
section IVIII 



B. Topological order 

To fully characterize our model, we must also ana- 
lyze the topological order of the fractionalizcd insulator 
H (|37p . We first recall the concept of topological or- 
der in 2D systems j^"i^^ Topologically ordered systems in 
two dimensions have three important physical properties. 
First, these systems have a finite energy gap separating 
the ground state(s) from excited states. Second, when a 
topologically ordered system is defined in a torus geome- 
try (periodic boundary conditions in both directions), the 
ground state is typically multiply degenerate. This de- 
generacy is not a consequence of any symmetry and is ro- 
bust to arbitrary perturbations)^^— Third, and perhaps 




FIG. 4. We can move a charge at position so around a flux 
at position Pq by applying a string of Uss' operators along a 
closed path C — sosi...SfeSo encircling Pq. 

most importantly, these systems have quasiparticle exci- 
tations with fractional statistics: when one such quasi- 
particle is braided around another, it acquires a nonvan- 
ishing Berry phase. 

Before analyzing the fractionalizcd insulator H (|37p . 
it is useful to first understand the topological order in 
the lattice boson model Hb ([I]). As we discussed earlier, 
this model contains two types of elementary excitations: 
"charge" excitations where gs = 1 at some site s, and 
"flux" excitations where Bp = e^'^'/™ on some plaquette 
P. We will now show that these excitations have non- 
trivial mutual statistics: when a charge moves around a 
flux, it acquires a Berry phase of 0ch,fi = 27r/m. 

Imagine we have a flux on plaquette Pq ^^nd a charge 
on some site sq. We can denote this state by 

|so,Po)-|n.„a = l,fepo = l) • (60) 

We need to adiabatically move the charge around the flux 
and compute the resulting Berry phase. To this end, we 
note that the operator Usqsi moves the charge from site 
So to a neighboring site si. 

Us.sAs^^Pq) = \si,Pq) (61) 

Indeed; this follows from the commutation relations be- 
tween Qt,Bp and Uss'- Therefore, we can move the 
charge around the flux by applying a string of Uss' oper- 
ators along some closed path C = soSi...SkSo encircling 
Po (see Fig. HI): 

Us,so-Us,s,Us„sAsO,Po) (62) 

To find the accumulated Berry phase, we use an operator 
identity: 

Us,so-Us,s,Usos^\so,Po) = n Bp\so.Po) 

Pec 

^e^^'/'"\so,Po) (63) 

Here, the second equality follows from the fact that bp ~ 
1 everywhere except for P ~ Pq, where bp^ = g27rj/m^ 

To isolate the statistical Berry phase from other geo- 
metric phases, we need to compare this phase with the 
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phase accumulated when the flux is not enclosed by the 
path C. In that case, the same operator identity gives 



(64) 



Comparing (|63| . (|64)) . we conclude that the statistical 
Berry phase is precisely 0ch,fl = 27r/m, as we claimed. 

In a similar way, we can check that there is no Berry 
phase associated with exchanging a pair of charge or flux 
excitationsi^ In other words, the charge and flux excita- 
tions are bosons. To complete our analysis of the topo- 
logical order, we need to find the ground state degeneracy 
on a torus. We describe this calculation in Appendix I A 21 
There, we show that the lattice boson model has exactly 
D = m? degenerate ground states in a torus geometry. 

The above ground state degeneracy and quasiparticle 
statistics are identical to the statistics and ground state 
degeneracy of the generalized "toric code" model with 
gauge group G = Zra^ Thus, the fractionalized bosonic 
insulator Rb (HI) has the same topological order as the 
Z„j toric code, or equivalently 2D gauge theory cou- 
pled to bosonic matter. 

Given these results, we can now easily analyze the 
topological order of the fractionalized insulator H ([57)) . 
This model contains three types of elementary excitations 
- charges, fluxes, and fermions. The charge and flux exci- 
tations are identical to the excitations in the lattice boson 
model and obey the same statistics. On the other hand 
the fermion excitation is a composite of an electron and 
k "charge" excitations, as we argued in equation p3p . 
This decomposition implies that the fermion excitation 
has mutual statistics 0f_fl = 2-Kk/m with respect to the 
flux excitations, but no mutual statistics with respect to 
the charges. As for the ground state degeneracy on a 
torus, it is easy to check that D = just as in the 
lattice boson model. 



C. Chern-Simons field theory description 

In this section, we derive a field theoretical description 
of the fractionalized insulator H (P7|) . This field theory 
is useful as it captures all of the (universal) properties of 
the fermionic phase-including both the topological order 
in the bulk, and the edge modes at the boundary (see 
section IVircTll . 

We begin by writing down a field theory description of 
the lattice boson model Hb H]). As we discussed in sec- 
tion [iVBl the topological order in this model is the same 
as Zm gauge theory coupled to bosonic matter. Previous 
work^ has shown that this kind of topological order is 
described by the Chern-Simons theory, 

Lb = ^e^^''{axdf,f3, + /3xdf,a,)-^e^^'''A^d^l3, (65) 
47r Ztt 

where Ax is the electromagnetic gauge field. (This par- 
ticular form of U{1) xU{l) Chern-Simons theory is some- 



times referred to as 2D BF theory^ii^^) . In this descrip- 
tion, the boson number current is given by 



— —AtJ-t^f) ft 

■Jboson o_ 

Zn 



(66) 



Also, the two types of quasiparticle excitations, "charges" 



{qs = 1) and "fluxes" {bp 



, are described by 



coupling (j65p to bosonic particles carrying unit and 
charge respectively. 

This field theory correctly describes the fractional 
charges and statistics of the quasiparticle excitations in 
the lattice boson model, as well as the ground state de- 
generacy on a torus. We can verify this using the "K- 
matrix" formalism2£i2ii2^ for abelian Chern-Simons the- 
ory and abelian fractional quantum Hall states. First we 
write Lb in ii'-matrix notation: 



Lb ^ " 



47r 



aixd^^aj, - ^tie^'^'Axd^ai, (67) 
Zn 



where 



Kij = 



m 
m , 



ti 



ai = 



(68) 



In this notation, the gauge charge carried by the two 
types of particles (charges and fluxes) can be represented 
in terms of the vectors 



(69) 



According to the ii'-matrix formalism, the physical elec- 
tric charge of each excitation is given by 



qi = l^K-H 



(70) 



while the mutual statistics associated with braiding one 
particle around another is given by 



(71) 



(The statistical phase associated with exchanging two 
identical particles is 6i = 9u/2). Applying this to Ich 
and /fl, we see that 

Qch = — , ffl = , 0,h.fi = — , 0ch = ^fi = (72) 
m m 

in agreement with the properties of the lattice boson 
model. Furthermore, according to the if-matrix formal- 
ism, the ground state degeneracy on a torus is given by 



D = \det{K)\ = 



(73) 



Again, this is in agreement with the properties of the 
lattice boson model. We conclude that the Chern-Simons 
theory ([S5)) does indeed describe the topological order in 
the lattice boson model. 
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To obtain a field theory description of tlie fractional- 
ized insulator H (|37p . we think of this system as a non- 
interacting topological insulator built out of fractionally 
charged fermions. These fermions are composite particles 
composed of k "charge" excitations of the lattice boson 
model together with one electron. In terms of the field 
theory (|65|) we can describe these composite particles by 
coupling (j65p to fermionic particles, each of which carries 
k units of charge and e units of electric charge. We 
conclude that the Lagrangian for H is of the form 



L ~ Lb + LTi[fca/, 



(74) 



where Lxiifco;^ + eA^J is the Lagrangian of a non- 
interacting topological insulator coupled to an external 
gauge field fca^ -)- eA^. 

To complete our derivation, we need to construct a 
field theory Ly/ for the non-interacting topological insu- 
lator coupled to an external gauge field. We can write 
down such a field theory explicitly in the case of the band 
structure (|49|) . In this case, the spin- up fermions form a 
band insulator with Chern number = +1, while the 
spin-down fermions form an insulator with Chern num- 
ber I' = — 1, so the appropriate field theory^! is a sum of 
two decoupled Chern-Simons theories: 



1 

'2^ 



(75) 



In this description, the spin-up and spin-down fermion 
currents are given by 



1 



(76) 
(77) 



while the spin-up and spin-down fermion excitations can 
be described by coupling (j75|) to bosonic sources carrying 
unit and 74^^ charge respectively. 

Combining (pSI) . (|74p . ([75)1 . we see that the fractional- 
ized electronic insulator H is described by the 4 compo- 
nent Chcrn-Simons theory 



L 

where 



47r 



aixd^aj, - — t/e^^'^AAa^a/, (78) 



2tt 



types of excitations, as well as the ground state degen- 
eracy on the torus. The "charge" excitations correspond 
to particles with unit ai charge, the "flux" excitations 
correspond to particles with unit 02 charge, and the spin- 
up/spin-down fermions correspond to particles with unit 
03 or 04 charge. We represent these particles with the 
vectors 







(0] 











1 


- h = 











1 






\p) 










V-V 



(80) 



Applying the formulas (|70ll7ip . we can check that the 
quasiparticle charges and statistics are 



'ch.fl 



2e 

m 
27r 



, <Zf ^ 
27rfc 



e(2fc 4- m) 



(81) 



7ch 



03 = 



in agreement with our previous calculations. Also, the 
ground state degeneracy on a torus is given by 



D = \det{K)\ 



(82) 



as expected. 



V. LATTICE MODELS FOR 3D FRACTIONAL 
TOPOLOGICAL INSULATORS 

In this section we construct exactly soluble lattice 
models for 3D fractionalized insulators. The 3D con- 
struction is a simple generalization of the 2D case, and 
much of our analysis carries over without any change. 
As before, we proceed in three steps, beginning by con- 
structing lattice boson models with fractional charge. 



A. Step 1: 3D lattice boson models with fractional 
charge 

The 3D lattice boson models are identical to the 2D 
models ^ , except that we consider these systems on the 
cubic lattice instead of the square lattice. In other words, 
the Hamiltonian is still 



Hs = vY^Ql~lY.^Bp + By) 



(83) 
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Kij = 



(79) 

As in the bosonic lattice model case (|67p . this Chern- 
Simons theory contains all the information about the 
fractional charge and fractional statistics of the three 



but now the sums over s and P run over the sites s and 
plaquettes P of the cubic lattice (Fig. Elja)). 

This model is exactly soluble just as in the 2D case, 
since Qs,Bp,Bp all commute with one another. Again, 
we can choose simultaneous eigenstates \qs,bp) where qs 
is an integer and 6p is a mth root of unity, and these 
states have energies 



u \ ^ 



-i,}ibp + b*p) 



2 ^ 

p 



(84) 
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FIG. 5. (a) The lattice boson model defined on the cubic 
lattice. As in the square lattice case, the Hamiltonian Hb 
(|83|l is a sum of Qs operators and Bp operators, but now the 
Qa operators act on six links (ss') and one site s, and the 
Bp operators act on the plaquettes of the cubic lattice, (b) 
The Bp operators obey the identity Hpgc ~ ^ where the 
product runs over the six plaquettes P adjacent to a "cube" 
C, and where we choose appropriate orientations on these 
plaquettes. 

Furthermore, as before, we can find the degeneracy of the 
Qs , bp eigenspaces for particular system geometries - say 
a rectangular slab of cubic lattice with open boundary 
conditions. The only difference from the 2D case is that 
the Bp operators satisfy the identity 

YlBp = l (85) 

Pec 

where the product runs over the six plaquettes P adja- 
cent to a "cube" C in the cubic lattice, and where we 
choose appropriate orientations on these plaquettes (Fig. 
[SJb)). Therefore, the bp quantum numbers obey the local 
constraint 

Hbp^l (86) 

for every "cube" C in the cubic lattice, in addition to the 
global constraint 

^ EE (mod m) (87) 

s 

We can check that for each set oi {qs,bp} satisfying these 
constraints, there is a unique eigenstate \qs, bp) (see Ap- 
pendix HH). 

As in the square lattice case, the ground state of the 
3D cubic lattice model (with open boundary conditions) 
is the unique state with ~ 0,bp ~ 1, and there are two 
types of elementary excitations. The first kind of excita- 
tion are "charges" where qs = 1 at some site s. These ex- 
citations are very similar to the charge excitations in the 
2D model. In particular, they carry the same fractional 
charge 2e/m, as can be verified using the arguments in 
section IIII A 21 On the other hand, the second kind of 
excitation is slightly different from the 2D case: instead 
of particle-like "flux" excitations, the 3D model has "flux 
loop" excitations where bp = e^'^'/™ along some closed 
loop in the dual cubic lattice. The reason that these ex- 
citations are loop-like in this case is the constraint ((86|) 
which ensures that the flux is divergence free and hence 
must form closed loops. 



B. Step 2: 3D lattice electron models with 
fractional charge 

The next step, as in the 2D case, is to modify the 
3D boson model ([55]) by including additional spin- 1/2 
electron degrees of freedom at each site s. The resulting 
3D electron model is identical to the square lattice case 
(PS)) . The Hamiltonian is still 

H, = vY,Ql-^Y.^Bp + Bl)-^iY,ns. (88) 

s P sa 

except now the sums over s, P run over sites and plaque- 
ttes of the cubic lattice. 

As in the 2D case, we can choose simultaneous eigen- 
states \qs,bp,nsa-) with corresponding energies 

^ = ^ E - ^ E(^^ + ^p) - ^ E (89) 

s P sa 

There is a unique state for each choice of g,,, 6p, Usa sat- 
isfying the two constraints 

6p = 1 , E '^■^ + ^" E ^^"'^ = ^ (™°^ ™) (^^) 

Pec s sa 

(assuming a geometry with open boundary conditions). 

Assuming < for simplicity, the ground state is 
the unique state with qa ~ 0,bp ~ l,nsa = 0. There 
are three types of elementary excitations, just as in the 
2D case: "charge" excitations with = 1, "flux loop" 
excitations with bp = e^'^'/™ along some loop in the dual 
lattice, and spin-1/2 "fermion" excitations where Jig^ ~ 
1. The spin-1/2 fermion excitations are similar to the 
fermion excitations in the 2D model, and in particular 
they carry the same fractional charge gf = e(l -I- 2k /m) 
(|34p . and transform under time reversal in the same way. 



C. Step 3: Building candidate 3D fractional 
topological insulators 

The final step is to modify the lattice electron model 
([55)1 so that the fractionally charged fermions form a 
non-interacting topological insulator. As in the 2D case, 
we accomplish this by adding a new term, i/hop to the 
Hamiltonian: 

H = H, + i/hop (91) 
where ffhop is defined just as in the 2D case ([38|) : 

i^hop = - E itss'aa'cl,^,CsaUt, + h.C.) (92) 

{ss') 

Again, we assume that tss'aa' u,V so that the band- 
width of the fermion excitations is small compared with 
the gap to the bosonic excitations. Just as in the 2D 
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case, we can show that the low energy physics of H is 
described by the non-interacting fermion model 



Hr.fr = 



E 

{ss') 



{tss'aa'd\,^,ds„ + h.c.) - fl^dl^dsa (93) 



We can then complete the construction, just as in the 
2D case: we simply choose the hopping amplitudes iss'o-o-' 
and chemical potential /x so that -ffeff is a non-interacting 
topological insulator. The low energy physics is then de- 
scribed by a topological insulator built out of fraction- 
ally charged fermions. There are many possible choices 
for tss'aa', but to bc concrete, we consider a 2-orbital 
model on the cubic lattice, introduced by Ref. [2^ The 4 
fermions on each site can be expressed as a single vector 
(dsi-j-, dsii, ds2\, ds2i)- In this basis the hopping matrix 
elements can be expressed as a tensor product of matri- 
ces Ti acting on the orbital indices (1, 2) and ai acting on 
the spin indices (t,i)- We take 



Ci —x.y,z 

+h.c.] + mSs-s',0 {tx ® (94) 

where a,a' £ {(l,t), (2,t), (1,1), (2,i)} specify both the 
spin and orbital indices. The first term is a spin-orbit 
coupling whose sign differs for the two orbitals; the sec- 
ond term is a nearest-neighbor hopping from one orbital 
to the other, and the third term is an on-site hopping 
between the two orbitals. Taking ^ = 0, together with 
t > and A > 0, this yields a topologically insulating 
band structure, provided that 2t < \m\ < 6t^. 



VI. PHYSICAL PROPERTIES OF THE 3D 
MODELS 



In the previous section, we introduced an electron 
Hamiltonian H (1911) whose low energy physics is exactly 
described by a non-interacting topological insulator of 
fractionally charged fermions. We now describe the phys- 
ical properties of the resulting fractionalized insulator. 



A. Surface states 

One of the most important properties of 3D non- 
interacting topological insulators is that they have gap- 
less surface states which are protected by time reversal 
symmetry and charge conservation)^ In this section, we 
discuss the analogue of these surface states for the frac- 
tionalized insulator realized by H 

To begin, we review the structure of the surface states 
of the non-interacting topological insulator with band 
structure ([94)) . To be concrete, we consider this model 
in an x Ly x slab of cubic lattice which is periodic 
in the x, y directions and has open boundary conditions 
at ^; = and z = L^- (This geometry is sometimes 



referred to as a "Corbino donut"). This system is trans- 
lationally invariant in the x, y directions, so we can still 
define a 2D band structure as a function of kx,ky. It 
has been shown that the surface band structure for the 
tight binding model ([M]) has two Dirac points — one for 
each surface. The surface states have a particularly sim- 
ple structure if the chemical potential is tuned near these 
Dirac points. In that case, the low energy surface states 
on each surface arc described by the Dirac Hamiltonian 



^surf 



cj, [o-a/3 • (*V + eA)] Ci3 - pLC^^Ca 



(95) 



where a,/3 =t, 4, a-nd cr = (ctj,, ctj^, (Tz). 

The surface states for the fractionalized insulator H 
(|9ip are almost identical to the non-interacting case, 
since the low energy physics of H is exactly described 
by the free fermion model ifeff- The only difference is 
that the fermions excitations carry fractional charge gf. 
Therefore, the low energy surface Hamiltonian for H near 
the Dirac points is given by exactly (|95p. except with e 
replaced by qi (and c replaced by d): 



HsuYi = vdl^ [(Tap ■ (iV + qfA)] dp - iid^^d-o 



(96) 



We leave the discussion of the stability of these surface 
modes to section IVlIl 



B. Surface quantum Hall effect 

The presence of a single Dirac cone in the surface mode 
spectrum has an interesting consequence if we break time 
reversal on this surface, for example by attaching to it a 
thin magnetic film. Let us first review what this does to a 
non-interacting topological insulator with band structure 
(1941) . Assuming that the magnetic film is only weakly 
coupled to topological insulator, we can model its effect 
by including a Zeeman term Hz = Bdl^a^pdp in the 
surface Dirac Hamiltonian (|95p . This term opens up a 
gap in the Dirac spectrum of size 2\B\. If the chemical 
potential lies within this gap, then the surface states are 
completely gapped and the system exhibits a half-integer 
surface Hall conductivity^!^ 



TI 



e 

2h 



(97) 



The total Hall conductivity on the two surfaces is e'^/h 
or depending on whether the magnetic field B has the 
same or opposite signs at z = and z = L^- 

The behavior of the fractionalized insulator H (PT|) is 
similar. As in the non-interacting case, the magnetic film 
introduces a Zeeman term, thereby gapping the spectrum 
of the surface Dirac Hamiltonian (|M)) . If the chemical 
potential lies within this gap, then the system exhibits a 
quantized surface Hall conductivity. The only difference 
is that the fermionic excitations carry charge qt so the 
Hall conductivity is a half-integer in units of qf /h: 



<lf_ 
2h 



(98) 
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The total Hall conductivity on both surfaces is q^/h 
or depending on whether the magnetic field B has the 
same or opposite signs at z = and z — L^. The first case 
is particularly interesting, since in this case the total Hall 
conductivity is fractional. Hence, if we view our geometry 
as a quasi-2D system, taking the thermodynamic limit 
Lx,Ly oo while keeping finite, then H realizes a 
fractional quantum Hall state. The properties of this 
state can be derived using the same approach as section 



C. Charge acquired by a magnetic monopole 

Another important property of non-interacting topo- 
logical insulators is their response in the presence of a 
magnetic monopole. Imagine threading a microscopic 
solenoid into a topological insulator which is so small that 
it fits within a single plaquette of the lattice. If we insert 
a full flux quantum $o = hc/e, we can effectively con- 
struct a monopole/anti- monopole pair at the two ends of 
the solenoid. An interesting property of non-interacting 
topological insulators is that this process causes a half- 
integer charge 

gxi = {n + l/2)e (99) 

to be localized at each end of the solenoid^^"— . In other 
words, a magnetic monopole (or anti-monopole) in a 
topological insulator carries half-integer charge. (Here 
the integer n depends on the microscopic form of the 
Hamiltonian near the ends of the solenoid. As such, it 
can be varied without encountering a phase transition in 
the bulk). 

We consider the analogous question for the fractional- 
ized insulators (|9T|) . In this case, we have to be careful 
because the Dirac quantization condition for monopoles 
is modified due to the presence of fractionally charged ex- 
citations. Rather than requiring monopoles to carry flux 
which is an integer multiple of $0j the Dirac argument 
implies that monopoles carry flux which is a multiple of 
(e/e*)<i>o where e* is the smallest charged quasiparticle 
excitation in the system. In our case, the minimal charge 
e* is given by 




if m is odd 
if m is even 



When TO is even, the bosonic charge excitations carry 
the minimal charge e*; when m is odd, an excitation 
with minimal charge e* can be constructed by forming a 
composite of an electron and (1— to) /2 charge excitations. 

If one violates the Dirac condition and inserts a 
monopole/anti- monopole pair whose flux is not a mul- 
tiple of (e/e*)$o, the result is that the monopole is 
joined by a physically observable Dirac string to the 
anti-monopole partner. One way to see this is to cal- 
culate the ground state energy in the presence of the 



monopole/anti- monopole pair using the 3D analogue of 
the eigenvalue spectrum (|56p. Examining the bp term, 
we can see that a Dirac string with flux j$o will have 
an energy cost ^ uL{l — cos(4j7r/m)) proportional to 
its length L unless 2j /m is an integer - that is, unless 
j is a multiple of e/e* ()100|) . (We cannot eliminate this 
energy cost by changing the value of bp along these pla- 
quettes due to the constraint Opgc^^'-) Intuitively, the 
point is that, unless this condition is satisfied, a charge e* 
quasiparticle excitation will acquire a nonvanishing Berry 
phase when it travels around the Dirac string. 

Hence, the natural quantity in the fractional case is the 
electric charge carried by a flux {e/e*)(j)Q monopole. To 
compute this charge, we again use the fact that the low 
energy physics of H ([5T|) is described by the free fermion 
model (|93p . Because the fermions carry charge qi, this 
system behaves exactly like a non-interacting topological 
insulator with a flux {q{/e*)(f>Q monopolc. In this case, we 
know that (n+ l/2)gf/e* fermions are localized near the 
monopole, as we discussed above. We conclude that the 
flux (e/e*)(f>o monopole in the fractionalized insulator H 
carries electric charge 

g = (n + l/2)qf/e* (101) 

On the other hand, we could presumably trap any 
number of additional charge e* quasiparticles near the 
monopole, by adding an appropriate local potential. 
Therefore, the only quantity which is independent of mi- 
croscopic details is the charge of the monopole modulo 
e* . This quantity is given by 

q=l^ (mod e*) if qt/e* is odd ^^^^^ 
1 (mod e*) if qi/e* is even 

More explicitly, if we use the expressions for qi p4p and 
e* poop , we find that qf/e* is even if and only if 2A: -f m 
is divisible by 4. 



D. Topological order 

In the previous two sections, we showed that the 3D 
model (PT|) exhibits a fractional surface Hall conductiv- 
ity (|M)) and a fractional charge bound to a magnetic 
monopole (|102p . These phenomena are manifestations 
of a fractional magnetoelectric effect^i^. On the other 
hand, it was argued in Ref. [l^ that the only way a frac- 
tional magnetoelectric effect can be consistent with time 
reversal symmetry is if the ground state on a three di- 
mensional torus is multiply degenerate. In this section, 
we show that our model does indeed have multiply de- 
generate ground states on a 3D torus. This ground state 
degeneracy originates from the fact that the model is 
topologically ordered, as we now discuss. 

Wc begin by analyzing the topological order of the 3D 
lattice boson model ([55]) . To derive the quasiparticle 
statistics in this model, we recall that (|83p contains two 
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types of elementary excitations: point-like "charge" ex- 
citations where = 1 at some site s, and "flux-loop" ex- 
citations where Bp = e^Tri/m along some loop in the dual 
lattice. Using the same approach as in the 2D case (sec- 
tion |IVB]), one can show that when one braids a charge 
around a flux-loop, it acquires a Berry phase of e^"^^^/™- 
where the sign depends on the orientation of the braiding 
trajectory. We can think of this as the 3D analogue of 
mutual statistics. It is also straightforward to show that 
there is no Berry phase associated with exchanging two 
charges, so the "charge" excitations are bosons. 

In addition to nontrivial statistics, the lattice boson 
model ((83|) also displays the second signature of topologi- 
cal order, namely ground state degeneracy. In particular, 
if we define this model in a 3D torus geometry, we find 
degenerate ground states (see Appendix IB 2\ . 

The quasiparticle statistics and ground state degener- 
acy of (j83|) is identical to the the statistics and ground 
state degeneracy of the 3D toric codoiii^ with gauge 
group G = Zjn- Thus, the lattice boson model (|83|) has 
the same topological order as the Z„i toric code in three 
dimensions, or equivalcntly 3D Z,„ gauge theory coupled 
to bosonic matter. 

Given these results, it is easy to analyze the topolog- 
ical order in the fractionalized insulator H ((9T|) . This 
model contains three types of elementary excitations - 
charges, fiux loops, and fermions. The charge and flux 
loop excitations are identical to the excitations in the lat- 
tice boson model and obey the same statistics. On the 
other hand the fcrmion excitation is a composite of an 
electron and k "charge" excitations, just as in the 2D 
case (j33p . This decomposition implies that the fermion 
excitation has mutual statistics e^'^*'^/™ with respect to 
the flux loop excitations. As for the ground state degen- 
eracy on a 3D torus, we can show that this is , just as 
for the lattice boson model. 



VII. GAPLESS BOUNDARY MODES: 
PROTECTED OR UNPROTECTED? 

In the previous sections we have constructed fractional- 
ized insulators (|37l9ip whose low energy physics is equiv- 
alent to a non-interacting topological insulator built out 
of fractionally charged fermions. We have also seen that 
these models have gaplcss boundary modes. However 
these two properties are not by themselves sufficient to 
declare these models "fractional topological insulators" : 
we must also show that the boundary modes are pro- 
tected. That is, we must show that the modes are robust 
to arbitrary time reversal invariant, charge conserving 
perturbations. 

In this section we investigate this robustness, and find 
two key results. First, we show that in both 2 and 3 
dimensions, the models for which the ratio qf/e* is odd 
have protected boundary modes. Thus, these models are 
indeed authentic "fractional topological insulators." Sec- 
ond ~ and perhaps more unexpectedly - we show that in 



the 2D case, the models with even qf/e* do not have 
protected edge modes, and are thus "fractional trivial in- 
sulators." We do not have an analogous result in the 3D 
case. That is, wc do not establish either the stability or 
instability of the surface modes when q^/e* is even. 

In the 2D case, our analysis of edge modes closely fol- 
lows the arguments of Ref. [l^. Our conclusion is also 
very similar to Ref. [loi : in that paper, the authors found 
that a certain class of conserving models have pro- 
tected edge modes if and only if asH/e* is odd, where 
UsH is the spin-Hall conductivity in units of e/27r, and e* 
is the elementary charge in units of e. Here, we find that 
our exactly soluble models have protected edge modes if 
and only if qi/e* is odd. The two criteria agree exactly, 
since as one can easily check, qf/e* = cTsh /e* (|50p for the 
models discussed here. 



A. 2D flux insertion argument 

In this section we present a general argument that 
shows that when qf/e* is odd, the 2D fractionalized in- 
sulator H ([57)1 has protected gapless edge modes. The 
argument is very similar to the one given in Ref. [Tol . 
which is in turn a generalization of the flux insertion ar- 
gument of Ref. H. 

To be precise, we prove the following statement: we 
consider the exactly soluble model 77 in a cylindrical ge- 
ometry with an even number of electrons and zero flux 
through the cylinder. We assume that the ground state 
is time reversal invariant and does not have a Kramers 
degeneracy on either of the two edges. (Wc discuss the 
meaning of this Kramers degeneracy assumption in sec- 
tion IVTlXT] below). Given these assumptions, we show 
that if qf/e* is odd, then there is always at least one ex- 
cited state whose energy gap vanishes in the thermody- 
namic limit. Furthermore, this excited state is robust if 
we add an arbitrary time reversal invariant, charge con- 
serving, local perturbation to the Hamiltonian, as long 
as we do not close the bulk gap. Finally, this excited 
state has the important property that it is in the same 
"topological sector" as the ground state, as we explain 
below. We interpret this low lying excited state as evi- 
dence for a protected gapless edge mode - that is, a mode 
which cannot be gapped out without breaking time rever- 
sal symmetry or charge conservation, explicitly or spon- 
taneously. 



1. Non-interacting case 

We flrst review the argument for the case of non- 
interacting topological insulators^ The proof that the 
system always has a low lying state begins as follows: 
we consider the system on a cylinder, and start in the 
many-body ground state \l/o at zero flux. We imagine 
adiabatically inserting $o/2 = hc/2e flux through the 
cylinder. Wc call the resulting state ^f^. Similarly, we 
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FIG. 6. The flux insertion argument in the 2D non- 
interacting case: we start in the ground state ^Jq and adi- 
abaticaUy insert ±"I>o/2 flux through the cyhnder, obtaining 
two states ^±tt. The state has a Kramers degeneracy at 
the two ends of the cyhnder, and is thus degenerate in energy 
with three other states, one of which is ^'-,r- If we start in 
one of these three degenerate states and then adiabatically 
reduce the flux to 0, we obtain an excited state \l'ex whose 
energy gap vanishes in the thermodynamic hmit. 



let Vf-TT be the state obtained by adiabatically inserting 
-$o/2 flux (See Fig. H]). We note that vj/^ and arc 
time reversed partners. 

To proceed further, we recall that 2D topological insu- 
lators are characterized by a Z2 invariant v = —l. This 
value of the invariant means that \l/7r has a Kramers de- 
generacy at the two ends of the cylinder if and only if 
^"0 does not.— (At an intuitive level, when we say that a 
many-body state has a Kramers degeneracy at the two 
ends of the cylinder, we mean that it can be brought 
to a state which is time reversal invariant by removing 
an odd number of electrons from each of the two ends. 
Assuming the ends are far apart, we can treat them as 
separate systems, each of which must have a two-fold de- 
generacy according to Kramers theorem. A precise def- 
inition of this notion of "local Kramers degeneracy" is 
given in Ref. [13) • By our assumption above, V^o has no 
Kramers degeneracy so ^'^ must have a Kramers degen- 
eracy at the two ends of the cylinder. Hence, as long as 
the two ends are well separated, Kramers theorem guar- 
antees that ^I^Tr is part of a multiplet of four states (two 
associated with each end) which are nearly degenerate in 
energy. More precisely, these four states are separated by 
an energy splitting which vanishes exponentially as the 
distance between the two ends grows. We note that vP-tt 
is one of these degenerate states, being the time reversed 
partner of 

Next, we imagine starting with the system at $o/2 
flux and then adiabatically decreasing the flux to 0. If 
we start with the state then adiabatic flux removal 
takes us to the ground state ^'o- However, if we start 
with Vf-T (or one of the other two degenerate states), the 
result is a eigenstate ^'ea; of the zero flux Hamiltonian, 
which is distinct from (see Fig. At the same time, 
it necessarily has low energy since the energy change AE 
associated with an adiabatic insertion of flux through 
a cylinder must vanish in the thermodynamic limit as 
long as charge conservation is not broken (in a gaplcss 



system with linear dispersion, we expect a scaling like 
AE 1/L, while in a gapped system, AE ~ (>-co7ist.-Ly 
To complete the argument, we now imagine adding an 
arbitrary time reversal invariant, charge conserving, local 
perturbation to the system (for example, wc could add 
short-ranged interactions between electrons). As long as 
the perturbation does not close the bulk gap, the above 
picture must stay the same: the Kramers degeneracy be- 
tween ^'tt and ^'-TT must remain intact, and hence ^'e^ 
must continue to be low in energy. We conclude that 
the system always contains at least one low-lying excited 
state whose energy gap vanishes in the thermody- 
namic limit. 



2. Fractionalized case 

We now explain the analogous argument for our frac- 
tionalized insulator H (|37p. The crucial difference from 
the non-interacting case is that this model has excita- 
tions with fractional charge e*. As a result, we need to 
modify the flux insertion argument, inserting ±^ • ^ 
flux through the cylinder instead of fluxiiS 

To understand why this is so, imagine we insert ±<f>o/2 
flux instead. Then the resulting states ^tt, ^'-tt differ by 
the insertion of a single flux quantum. But the inser- 
tion of a single flux quantum changes the Berry phase 
associated with moving a charge e* particle around the 
cylinder by Ad = 2tt^. As Ad is not a multiple of 27r, 
^Ptt and Vt-TT are in different "topological sectors." That 
is, to get from one state to the other, we need to transfer 
a fractionalized excitation (in this case, two flux quasi- 
particles - see section UlIPp from one end of the cylinder 
to the other jiSii^ If we then try to construct a low ly- 
ing excited state at zero flux in the usual way (i.e. by 
starting from ^'_7r and adiabatically reducing the flux to 
zero), the resulting state ^'ex is in a different topologi- 
cal sector from ^fQ. Unfortunately, this property means 
that we cannot conclude much from the existence of ^'e^. 
The reason is that topologically ordered systems like H 
always have a finite number of low lying states that lie in 
different topological sectors.—"— These states have noth- 
ing to do with time reversal symmetry and continue to 
exist even if time reversal symmetry is broken and the 
edge is gapped. Instead, these states are "topologically 
protected" and are analogous to the degenerate ground 
states in a torus geometry. Therefore, it is important 
that we construct a low-lying state in the same topolog- 
ical sector as the ground state. This will establish the 
existence of an "unexpected" low lying state, which can 
plausibly be taken as evidence for a time reversal pro- 
tected gaplcss edge mode. 

For this reason, wc insert zL^^-j^ flux through the cylin- 
der where TV = e/e*: when we insert this amount of flux, 
the states '^nttj'^-Ntt lie in the same topological sector, 
so we do not run into the issue discussed above. Other 
than this simple modification, the argument proceeds as 
before: wc note that the low energy physics of H is de- 
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scribed by a topological insulator built out of charge qi 
fermions, so our insertion of ±^ • *I'o/2 flux is equivalent 
to inserting ±p- • flux for a non-interacting topological 
insulator. If qf/e* is odd, then this process changes the 
time reversal polarization at the two ends of the cylin- 
der, implying that ^n-r^'^'-Ntt are Kramers degenerate. 
As before, we can then construct a protected low lying 
state at zero flux by starting with ^'-at^ and adia- 
batically inserting flux quanta. Furthermore, ^ea; 

is in the same topological sector as the ground state V^o 
since '^'ntt, '^-Ntt are in the same topological sector. We 
conclude that when qi/e* is odd, the fractionalized in- 
sulator H has a protected low lying excited state in the 
same topological sector as the ground state. This is what 
we wanted to show. 



B. 3D flux insertion argument 

Next, we show that a similar flux insertion argument 
can be used to establish the existence of protected low- 
lying surface states in the 3D models (PT|) with qi/e* odd. 
The precise statement we will prove is this: we consider 
the 3D models in a "Corbino donut" geometry, periodic 
in the x and y directions, and with open boundary condi- 
tions in the z direction (see Fig. [Tji) . Wc assume that the 
number of electrons is even and that the ground state is 
time reversal invariant when there is zero flux through the 
two holes of the Corbino donut. In particular, we assume 
that the ground state does not have a Kramers degener- 
acy on the z = 0, z ~ surfaces. We then show that 
if qt/e* is odd, then there is always at least one excited 
state whose energy gap vanishes in the thermodynamic 
limit. Furthermore, this excited state is robust if wc add 
an arbitrary time reversal invariant, charge conserving, 
local perturbation to the Hamiltonian, as long as we do 
not close the bulk gap. Finally, this excited state is in 
the same topological sector as the ground state. Just as 
in the 2D case, we interpret this low lying excited state 
as evidence for a protected gapless surface mode-that 
is, a mode which cannot be gapped out without breaking 
time reversal symmetry or charge conservation, explicitly 
or spontaneously. 



1. N on-interacting case 

It is useful to flrst give the argument in the case of the 
non-interacting 3D topological insulator. (This argument 
was alluded to in Ref . though not explicitly discussed) . 
The proof that the system has a protected low lying state 
begins similarly to the 2D case: we start in the many- 
body ground state at zero flux, which we call *I'(o.o)i and 
then imagine adiabatically inserting $0/2 flux through 
each of the two holes of the Corbino donut. In this way, 
we obtain three new states ^'(7r,o)j ^(o,7r); ^(tt.tt)- 

We next recall that 3D topological insulators arc char- 
acterized by a Z2 invariant v = —1. Similar to the situa- 
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FIG. 7. (a) The 3D flux insertion argument makes use of a 
"Corbino donut" geometry, with flux <j)x,(l>y through tlie two 
lioles of tlie Corbino donut. (b) In the non-interacting case, 
we start in the ground state ^^(o.o) a-nd adiabatically insert 
$0/2 flux through each of the two holes, obtaining 3 other 
states *I'(,r,o)) *I'(o,7r)5 'I'(ir,7r)- Au odd number of these states 
must have Kramers degeneracies at the two surfaces of the 
Corbino donut, leading to four possible degeneracy patterns 
(filled circles represent states with Kramers degeneracies). In 
particular, at least one state has a Kramers degeneracy, and 
hence can be used to construct a protected low lying excited 
state "^ex, as in the 2D case. 



tion in 2D, this value of the invariant implies that an odd 
number of the 4 states ^'(o,o)j 'J'(7r,o)j ^(o,7r)j ^(7r,7r) have 
a Kramers degeneracy on the z = and z = surfaces 
of the Corbino donut 1^ (We give a simple derivation of 
this result at the end of this section). By our assump- 
tion above, ^'(o.o) has no Kramers degeneracy, so either 
all three of the other states have a degeneracy, or ex- 
actly one of them does (see Fig. Wp)- In particular, at 
least one of the three other states has a Kramers degen- 
eracy. Let us denote this state by '^(a,p)- Then, as long 
as Lz is large so that the two surfaces of the Corbino 
donut are well separated, Kramers theorem guarantees 
that (a,i3) is part of a multiplct of four states (two as- 
sociated with each surface) which are nearly degenerate 
in energy. More precisely, these four states are separated 
by an energy splitting that vanishes exponentially as the 
distance between the two surfaces grows. 

Just as in the 2D case, we can use this Kramers degen- 
erate state to construct a low lying excited state at zero 
flux: as in that case, wc start in one of the three states 
which are degenerate with ^[a^p), and then adiabatically 
reduce the flux to (0, 0). The result is an eigenstatc ex 
of the zero flux Hamiltonian. This state is distinct from 
the ground state ^'(o,o) and its energy gap vanishes in 
the thermodynamic limit. Furthermore, time reversal in- 
variant perturbations cannot change this picture, since 
they cannot split the degeneracy between ^'(q./j) and its 
Kramers partners. Hence, even in the presence of arbi- 
trary time reversal invariant, charge conserving pertur- 
bations, the system always has at least one low lying 
excited state '^ex- 

For completeness we now explain why non-interacting 
topological insulators have the property that an odd 
number of the 4 states ^(o,o)j ^(ir,o)j ^(o,ir)j ^(tt.tt) have 
a Kramers degeneracy on the z = and z = surfaces 
of the Corbino donut. This result was first established 
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then depending on the parity of Lx,Ly, the momentum 
states near the Dirac point at flux {(f>x-,4'y) = (0,0), can 
be arranged in any of the three possibihties shown in Fig. 
IHt-d. Again, ^"(0.0) will be a Slater determinant involving 
all states in the lower band. In these three cases, one 
can check that when one inserts flux through the two 
holes of the Corbino donut, exactly one of the resulting 
states 4'(7r,o)i ^(o,7r)5 ^(ir.Tr) has an unpaired momentum 
state, and therefore a Kramers degeneracy on the = 
surface. Combining all of these cases, we conclude that 
the system always has a Kramers degeneracy in either 1 
or 3 of the four possible states. This is what we wanted 
to show. 



FIG. 



Surface state configurations 



of 



'I'(o,o),*(7r,o),*(o,7r),'I'(,r,,r) for a systsm with a single 
Dirac cone on each surface, and a chemical potential just 
above the Dirac point. Here we plot the surface states as a 
function of kx', the bands indicate different discrete values 
of ky. The band that crosses the tip of the Dirac cone, 
corresponding to ky — 0, is singly degenerate; all other 
bands, for which ky 7^ 0, are doubly degenerate. The 4 
panels show the arrangement of the discrete surface momenta 
relative to the tip of the Dirac cone in each of the 4 flux 
sectors. The filled circles indicate occupied states. 



in Ref. |^; here we give an alternative derivation using 
the fact that these systems have an odd number of Dirac 
cones on each surface. 

For simplicity we assume a band structure with a single 
Dirac cone on each surface. We also assume that the 
chemical potential is tuned to lie just above the Dirac 
point. We begin with the system at zero flux: {(j)x, (py) = 
(0,0). Letting and Ly be the dimensions in the two 
periodic directions, the momenta kx, ky are quantized as: 



2'Kmx 



2'Kmy 

Ly 



(103) 



where nixTmy are arbitrary integers. 

Let us focus on one of the two surfaces, say the z — 
surface. Suppose that the Dirac cone on that surface is 
located at {kx,ky) = (0,0). Then the momentum states 
near the Dirac point will be arranged as in Fig. [5^. The 
wave function ^'(o,o) "^iH be a Slater determinant consist- 
ing of all the momentum states with < /i - including 
both momentum states at {kx,ky) = (0,0). If we now 
thread half a flux quantum through the x or y direc- 
tion, the momenta will shift by half a unit in kx or ky, 
leading to the states vE'(7r,o): ^(o,7r), ^'(•n-.Tr) shown in Fig. 
Hh-d. Examining these configurations, we can see that all 
three have a single filled state in the upper band, imply- 
ing that all three of them have a Kramers degeneracy on 
the 2; = surface (these states must also have a Kramers 
degeneracy on the z = surface since the total number 
of electrons is even). 

On the other hand, if the Dirac cone is located at a 
different point in the Brillouin zone, say (kx, ky) = (tt, tt). 



2. Strong vs. weak topological insulators 

The reader may have noticed that, in the flux inser- 
tion argument, we only used the fact that at least one of 
the states ^^(Tr.o): ^(o,ir)j 'J'(7r,7r) had Kramers degenera- 
cies: we did not make use of the additional information 
that the number of such states was odd. In other words, 
the argument would have worked equally well if 2 of these 
states had degeneracies. 

We can understand the physical meaning of this obser- 
vation as follows. Recall that, in addition to the usual 
3D topological insulator, there is another interesting kind 
of 3D time reversal invariant band insulator known as a 
"weak topological insulator."— (In this terminology, the 
usual 3D topological insulator is known as a "strong topo- 
logical insulator.") Weak topological insulators can be 
thought of as layered systems, where each layer is a 2D 
topological insulator. The reason that weak topological 
insulators are relevant here is that these systems give ex- 
amples where 2 of the three states ^'(,r,o)i ^(o,ir)j ^(ir.;r) 
have Kramers degeneracies (assuming, as before, that 
^(0.0) is not Kramers degenerate). More precisely, this 
pattern of degeneracies must arise in a weak topological 
insulator with an odd number of layers, with the layers 
oriented perpendicular to the z ^ 0,Lz surfaces. 

This example clarifies why the fiux insertion ar- 
gument extends to systems where 2 of the states 
^(tt.o), *(o,7r)> ^(7r,7r) have Kramcrs degeneracies: just like 
strong topological insulators, these odd-layer weak topo- 
logical insulators have protected surface modes. One way 
to understand these surface modes is to note that each 
2D topological insulator layer contributes a one dimen- 
sional edge mode at the boundary. These edge modes can 
be gapped out in pairs by appropriate perturbations, but 
if there are an odd number of layers, we will always be 
left with at least one gapless mode. 

At the same time, this example reveals a shortcom- 
ing of the flux insertion argument: the argument does 
not distinguish strong topological insulators from weak 
topological insulators with an odd number of layers, as 
both systems have protected surface modes. Yet it is 
clear intuitively that the surfaces of these two systems 
are topologically distinct: the surface modes of a strong 
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topological insulator are in some sense "more two dimen- 
sional" than those of a weak topological insulator with 
an odd number of layers. 

In the following, we refine the flux insertion ar- 
gument so that it addresses this issue. We use 
the additional information that an odd number of 
^(0,0) J 'J'(ir,o)j ^(o,ir)j ^(tt.tt) have Kramers degeneracies to 
show that our system not only has a low lying surface 
mode, but this mode is delocalized in both the x and y 
directions. To be precise, we prove that the low lying 
surface modes in our system can never be localized to 
one dimensional strips. This property distinguishes our 
system from a weak topological insulator with an odd 
number of layers since in that case a carefully designed 
perturbation can gap out all low lying excitations ex- 
cept a single one dimensional edge mode. (As an aside, 
we note that the difference between the surface modes of 
weak and strong topological insulators largely disappears 
if one considers the surfaces in the presence of a random 
potential^). 

We give a proof by contradiction. Suppose that an 
appropriate perturbation could localize all low-lying sur- 
face modes to a single one dimensional strip. Let Wx, Wy 
denote the number of times this strip winds around the 
(toroidal) surface in the x and y directions respectively. 
Then, when we insert flux (j)x,4'y = 0, tt through the 
two holes of the Corbino donut, the low lying excita- 
tions will only couple to 4>x , (t>y via the linear combination 
0strip = Wx4>x + Wy(j)y In particular, it follows that the 
presence or absence of a Kramers degeneracy at '^(if,^^^^) 
only depends on (t>x,4'y via ^strip = Wx4>x + Wy(j)y. But 
this kind of 4>x , ipy dependence always leads to an even 
number of 1'(o,o), *(7r,o), ^(o.tt), ^CTr.Tr) having Kramers 
degeneracies, as we now show. 

It is useful to consider the cases where Wx , Wy are 
even and odd separately. First suppose that Wx is even 
and Wy is arbitrary. In this case, the effective flux 
0strip takes the same value at {4>x, (f'y) = (0, 0), (tt, 0) and 
similarly at {(f)x,(t>y) = (0, tt), (tt, tt). Therefore, since 
0strip completely determines whether the system has a 
Kramers degeneracy, it follows that either ^1/(0,0); ^(ir,o) 
both have Kramers degeneracies or they both have no 
degeneracy, and similarly for Vl'(o,7i-)j ^(ir,7r)- But then 
the Kramers degenerate states come in pairs, implying 
that an even number of the four states have Kramers 
degeneracies. The same is true if Wy is even and Wx 
is arbitrary. The only remaining possibility is if Wx , Wy 
are both odd. In this case, ^strip takes the same value 
at {(j}x,(j)y) — (0, 0), (tt, tt) and similarly at {(j)x,4>y) ~ 
(0, tt), (tt, 0). Again, the Kramers degeneracies come in 
pairs and we are guaranteed to get an even number of 
Kramers degenerate states. We conclude that a system 
where an odd number of ^'(0,0)7 '^{7^,0}^ ^(o,7r)i ^(tt.tt) have 
Kramers degeneracies must have excitations which are 
truly two dimensional and are not localized to a one di- 
mensional strip. 



3. Fractionalized case 

The argument for the 3D fractionalized insulator H 
(|9T|) is completely analogous. For the same reason as 
in 2D, we insert flux rather than $o/2 flux, 

where N = e/e*. Inserting flux through each of the 
two holes in the Corbino donut, we construct four states, 
*(o,o), '^{N7T,o), *(o,JVir), *(Af7r,Af7r)- Siuce the low energy 
physics is described by a topological insulator built out 
of charge qt fermions, this flux insertion is equivalent 
to inserting • ^ flux for a non-interacting topologi- 
cal insulator. If qf/e* is odd, then an odd number of 
*(o,o),*(Ar7r,o),*(o,JV7r),*(Ar7r,Ar7r) havc Kramers degen- 
eracies. We can then construct a protected low lying 
state ^'ea; at zero flux, by starting in one of the Kramers 
degenerate states and adiabatically reducing the flux to 
0. Moreover, this state is in the same topological sector 
as the ground state, just as in the 2D case. This is what 
we wanted to show. 

Finally, we note that we can refine the flux insertion ar- 
gument just as in the non-interacting case. To be specific, 
by repeating the analysis in the previous section, we can 
prove that if qf/e* is odd, the fractionalized insulators not 
only have low lying surface modes, but these modes are 
truly two dimensional and cannot be localized to a one 
dimensional strip. Just as in the non-interacting case, 
this allows us to distinguish these systems from "weak 
fractional topological insulators" with an odd number of 
layers. 



C. Microscopic theory of 2D edge 

While the flux insertion argument shows that the edge 
modes are protected for the models ([57]) with odd qi/e*, it 
does not tell us anything about the converse statement — 
that is, whether the modes can be gapped out in the 
models with even qf/e* . In order to establish this fact we 
now analyze the microscopic theory of the edge. 



1. Constructing the 2D edge theory 

The first step is to construct the edge theory. The free 
fermion theory ()59|) is unfortunately not a good starting 
point since it is not truly a complete edge theory. The 
reason it is not complete is that it doesn't include all the 
quasiparticle excitations with nontrivial statistics, such 
as flux quasiparticles. As a result, this theory cannot be 
used to analyze the most general edge perturbations — 
for example, those that are large compared to the flux 
quasiparticle gap. 

A simple way to construct a complete edge theory is 
to use the general bulk-edge correspondence for abelian 
Chern-Simons theory.— 1^ According to this correspon- 
dence, the edge of H ([57)) is described by a four compo- 
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ncnt chiral boson theory 



L 



1 
1 

2^ 



(104) 



Here $ is a four component vector of fields, V is the ve- 
locity matrix and K and t are defined as in ((79)) . In this 
language, the most general quasiparticle creation opera- 
tors can be written as e*' *, where / an integer valued 
four component vector. The charge carried by e*' * is 
l^K~^t. The subset of these operators that are "local" 
(i.e. products of electron creation and annihilation op- 
erators) can be written as e*®*^^^ where 9(A) = A^K^ 
and A is an integer valued four component vector. 

The final component of the edge theory is the phys- 
ical interpretation of the e'®'-^-' operators in terms of 
the microscopic model H. Using our understanding of 
the corresponding bulk Chern-Simons theory, we iden- 
tify A"^ = (0, 0, 1, 0), (0, 0, 0, 1) with creation operators 
for spin-up and spin-down electrons, respectively, and 
A-'" = (0, 1,0,0) with a creation operator for the charge 
2e spinless boson. As for (1, 0, 0, 0), this creates a neutral 
spinless particle which can be thought of as composite of 
m flux quasiparticles. The most general e*®'-^-' is a com- 
posite of these elementary operators. 

These identifications have the added benefit of fixing 
the transformation properties of $ under time reversal. If 
we require that the electron creation operators transform 
as 
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(105) 



and that the spinless charge 2e boson is invariant under 
time reversal, while the flux changes sign, we deduce that 
$ transforms as 



$ ^ T<I> + ttR-^x 



(106) 
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Equivalently, 8(A) transforms as 

e(A) ^ e(-rA) - 



(107) 



Q(A) 



(108) 



where Q(A) = A-^x- To complete the story, we need 
to explain the relationship between the above edge the- 
ory (|104|) and the edge ([59| of the exactly soluble model 
H ([57)) . An important clue is that ()104p has four gap- 
less modes while the free fermion edge ([59]) has only two 
modes. This mode counting suggests that ([59| can be ob- 
tained from ()104[) by adding a perturbation that gaps out 
two of the modes. In appendix [Cl we confirm this guess. 



We show that a time reversal invariant perturbation of 
the form 



U{A) = C/(a;)[cos(e(A) ~ a{x)) 

+ (-l)Q(A) cos{Q{TA) ~ aix))] 



(109) 



with Aq = (1,0,0,0) does the job. That is, the pertur- 
bation U{Ao) gaps out two of the edge modes of ()104p 
leaving behind exactly the free fermion edge ([59]). 

Turning this statement around, if we start with the 
exactly soluble model ([37| . we must be able to find a 
perturbation that closes the gap to the bosonic excita- 
tions at the edge, and results in the four component edge 
theory ()104|) . At a microscopic level, the following per- 
turbation accomplishes this task: 



AH, 



edge 



J E (Us 

{ss')edX 



h.c) 



(110) 



Here, the sum runs over links on the boundary dX of 
the exactly soluble system. Physically, this term gives 
an amplitude for the bosonic charge excitations to hop. 
When J <^V , the bosonic charge excitations arc gapped 
in both the bulk and the edge, and the only low energy 
edge excitations are the two chiral fermion modes in ()59|) . 
However, when J is sufficiently large, the gap to the 
bosonic charge excitations closes at the edge, resulting 
in a ID superfluid. The resulting edge has four gapless 
modes and is described by ()104p . 



2. Analysis of edge mode stability 



In the previous section we derived an edge theory ()104|) 
for H ([37]) , and showed it could be obtained from the ex- 
actly soluble edge by adding an appropriate perturbation. 
We now investigate whether this edge theory (|104p can be 
gapped out completely by charge conserving, time rever- 
sal symmetric perturbations. We find that the edge can 
be gapped out if and only if qi/ e* is even, in agreement 
with the flux insertion argument. 

Our analysis closely follows that of Ref. [T^. We focus 
on scattering terms of the form ()109p . These perturba- 
tions can be divided into two different classes: perturba- 
tions where A,rA are linearly independent, and pertur- 
bations where TA = ±A. Perturbations of the first type 
can gap out four edge modes, while perturbations of the 
second type can gap out two edge modes. Therefore, in 
order to fully gap out the edge ([104p . we either need one 
perturbation of the first type or two perturbations of the 
second type. Here, we will focus on the second possibility 
(though we would obtain the same results if we consid- 
ered the first kind of perturbation instead). That is, we 
will look for perturbations of the form 



C/(Ai) + C/(A2) 



(111) 



where TAi = ±Ai and similarly for A2. 

To begin, we note that the most general charge con- 
serving solution to TA = A is A-^ = (0, x, — a;, — x). 
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while the most general solution to TA = — A is A'^ = 
(y, 0, z, — z). Therefore, in order to get two linearly inde- 
pendent A's, we can either take 

Af = (yi,0,zi,-zi) , A^ = (y2, 0,^2,-2:2) (112) 

or we can take 

Af = (0, X, -X, -x) , = {y, 0, z, -z) (113) 

In the first case (|112p . the corresponding perturbation 
(jllip can indeed gap out the edge, but hand in hand with 
that it spontaneously breaks time reversal symmetry, as 
we now explain. We note that when (|llip gaps out the 
edge, it freezes the value of 0(Ai),0(A2). It therefore 
also freezes the value of 0(?/iA2 — y2Ai), which is a mul- 
tiple of 6(0, 0, 1, -1). But cos(e(0, 0, 1, -1) - a) is odd 
under time reversal (this follows from (|108l) or alterna- 
tively from the fact that cos(0(0, 0, 1, —1) — a) describes 
a Zeeman field oriented in the xy plane) so we conclude 
that this perturbation spontaneously breaks time rever- 
sal symmetry. 

Therefore, if we want to gap out the edge without 
breaking time reversal symmetry, we are led to pertur- 
bations of the form (|113|) . According to Haldane's null 
vector criterion22 such a perturbation can gap out the 
edge if Ai, A2 satisfy 

AJkAi = A^KA2 = AJKA2 = (114) 

(The origin of the criterion (|114p is that this condition 
guarantees that we can make a linear change of variables 
from $ to such that (a) the action for $' consists of 
two decoupled non-chiral Luttingcr liquids, and (b) the 
two perturbations correspond to backscattering terms for 
the two liquids. See appendix [Q for a related example). 
For the above Ai, A2, this leads to the condition 

(2fc + m)?;-2z = (115) 

It is convenient to parameterize y, z by y = us, z = vs 
where u, v have no common factors. Then the above 
condition becomes 

(2fc + m)u-2v = (116) 

We now consider two cases: 2k + m divisible by 4, and 
2k + m not divisible by 4. If 2fc -I- m is not divisible by 
4 then we must have v odd (since even v implies even 
u, which contradicts the fact that u,v have no common 
factors). It then follows that the perturbation U{A2) 
spontaneously breaks time reversal symmetry: the per- 
turbation U{A2) freezes the values of Q{A2) which is a 
multiple of Q{u, 0, v, —v). But cos(9(u, 0, v, —v) — a) is 
odd under time reversal symmetry (since v is odd), so 
this perturbation must spontaneously break time rever- 
sal symmetry. We conclude that when 2k + m is not 
divisible by 4, we cannot gap out the edge using these 
perturbations without breaking time reversal symmetry. 
We note that this agrees with the flux insertion argument 



since qt/e* is odd precisely when 2k + m is not divisible 
by 4. 

On the other hand, when 2fc + m is a multiple of 4 (or 
equivalently, qt/e* is even) the above condition suggests 
a natural solution for (Ai, A2): we can take u = 1, v = 

k + m/2, s = 1, a; = 1, so that 

Af = (0,1,-^1,-1) 

A^ = (l,0,fc-f m/2,-fc-m/2) (117) 

The physical meaning of the corresponding perturbations 
is clear: U{Ai) describes a process where a charge 2e 
boson breaks into two electrons with opposite spin di- 
rections, while U{A2) describes a process which flips the 
spins of fc 4- to/2 electrons while simultaneously creating 
a composite of to flux quasiparticles. 

To see that the perturbation U{Ai) + U{A2) can gap 
out the edge, we note that the above Ai,A2 obey Hal- 
dane's criterion (|114p . Moreover, this perturbation is 
time reversal invariant as long as ai = 7r/2 or 37r/2. 
Finally, this perturbation does not break time reversal 
symmetry spontaneously, as we now explain. The ba- 
sic point is that, the only way that time reversal sym- 
metry (or any other symmetry) can be spontaneously 
broken is if, for some a, b with no common factors, the 
linear combination aAi + 6A2, is non-primitive - that is, 
aAi + 6A2 = nA, where A is an integer vector, and n 
is an integer larger than 1. (To understand where this 
condition comes from, see the example of spontaneous 
symmetry breaking given in the discussion after (|113p ). 
But it is clear from the form of Ai , A2 that all such linear 
combinations oAi -I- 6A2 are primitive. We conclude that 
U (Ai) -I- U{A2) gaps out the edge without breaking time 
reversal symmetry explicitly or spontaneously. 

VIII. CONCLUSION 

In this work, we have constructed a family of exactly 
soluble models for fractionalized, time reversal invariant 
insulators in 2 and 3 dimensions. At low energies, these 
models behave like topological insulators made of frac- 
tionally charged fermions of charge qf. As a result, the 
2D models have two edge modes of opposite chiralities, 
while the 3D models have a gapless Dirac cone on each 
surface. At high energies, the insulators possess other 
types of excitations, including quasiparticles with mini- 
mal charge e*. As expected for systems with fractional 
charge, all of the 2D and 3D insulators are topologically 
ordered, exhibiting fractional statistics as well as ground 
state degeneracy in geometries with periodic boundary 
conditions. 

An important characteristic of the 3D models is that 
they exhibit a fractional magnetoelectric effect. More 
specifically, if time reversal symmetry is broken at the 
surface, the 3D models exhibit a surface Hall effect with 
fractional Hall conductivity a^y = qf/'^h. In addition, if 
a magnetic monopole is inserted into the bulk, it binds 
a fractional electric charge. Somewhat surprisingly, the 
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size of the elementary charge e* plays a significant role 
in this charge binding physics. 

A key question is whether the gaplcss boundary modes 
are robust to perturbations, so that these systems truly 
qualify as "fractional topological insulators." We have 
studied this question in detail and have shown that the 
boundary modes are protected for the models where the 
ratio Qi/e* is odd. That is, the edge or surface modes can- 
not be gapped out by perturbations that are time reversal 
symmetric and charge conserving. In fact, these modes 
are immune even to large perturbations, as long as the 
perturbations do not close the bulk gap or spontaneously 
break time reversal or charge conservation symmetry at 
the boundary. 

In contrast, the gapless boundary modes are not nec- 
essarily stable when gf/e* is even. In the 2D case, we 
have demonstrated this point explicitly by constructing 
time reversal invariant, charge conserving perturbations 
that open a gap at the edge. The situation in 3D is less 
clear: though our argument for proving the existence of 
protected surface modes breaks down in the models with 
even (jf/e*, we have not explicitly constructed perturba- 
tions that gap out the surface. Determining the robust- 
ness of the surface modes in these models is an interesting 
question for future research. 

Our findings regarding the stability of the edge modes 
arc intriguing for several reasons. First, we have demon- 
strated that simply having fractionally charged fermions 
in the right band structure does not guarantee a system 
is a fractional topological insulator. Second, although e* 
is a high energy property, it evidently determines the fate 
of the low energy excitations on the edge. And third, it 
is interesting that the condition that guarantees stability 
the boundary modes turns out to be identical for both 
the two and three dimensional models considered here 
and those considered in Ref. [l^. 

We have focused here on a particular set of models, 
which does not exhaust all the possibilities for fractional 
topological insulators. In two dimensions, the fractional 
quantum spin Hall systems discussed by Refs. andfiol 
give a whole other class of examples. Similarly, in three 
dimensions, there may be other families of 3D fractional 
topological insulators beyond the ones discussed here. 
Constructing microscopic models of the other families, 
and finding a complete classification of the possibilities, 
remain tantalizing open questions. 
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Appendix A: Eigenstate degeneracy of 2D lattice 
boson model 



In this section, we find the degeneracy of the \qs,bp) 
eigenstates of the 2D lattice boson model ([T|). We con- 
sider two geometries: a rectangular piece of square lattice 
with open boundary conditions (Fig. [H^a)), and a rect- 
angular piece of square lattice with periodic boundary 
conditions — that is, a torus (Fig. [2Ib)). 



1. Open boundary condition geometry 

The first step is to define projectors Pq^ , P^p , which 
project onto states with Qs = Qs and Bp ~ bp respec- 
tively. The degeneracy D of the qs,bp eigenspace can 
then be written as the trace of the product of all the 
projectors: 



D = Tr (j{P,.\{Pb)j 



(Al) 



We next write down an explicit expression for Y\p P^p . 
To this end, we recall that the eigenvalues oi Bp are pth 
roots of unity, so the projector can be written as 



-. m — 1 



(A2) 



3=0 

Expanding out the product, we have 



P P \ j=0 

-;^En"^i^^n^? (A3) 

where A'^piaq is the total number of plaquettes. Using the 
definition Bp = Ul2U2^Uz^Ul^l^ we get an expression of 
the form 



(A4) 



{jp} P 



{SS-) 



Here, Ajss' = jp—jp' where P, P' are the two plaquettes 
bordering the link {ss'). If {ss') happens to live on the 
boundary of the system, then Ajss' = jp where P is the 
unique plaquette bordering this link. Substituting this 
into (|Aip gives 



{jp} P \ {ss') J 

The next step is to note that the operator [/J^, changes 
the boson number Uss' by ±r (mod m). We conclude 

that the trace Tr Pg^ H{ss') ^.ss''""') vanishes unless 
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Ajss' = (mod m) everywhere, which implies jp = 
everywhere (since we are assuming a geometry with open 
boundary conditions). In this way, we see that the only 
nonvanishing term in (jASP is the one with jp = 0, so that 



D 



(A6) 



All that remains is to compute this trace, or equivalently, 
count the number of states with Qg = Qs- Working in the 
number basis, \ns,nss'), and using the definition ^ of 
Qs it is clear that there is a one-to-one correspondence 
between states with Qs = qs and configurations of Uss' 
satisfying 



qs (mod m) 



(A7) 



We can count these configurations by comparing the 
number of constraints (jA7p to the number of free param- 
eters. Examining (|A7p . we see that we have a constraint 
for every site s of the lattice. However, these constraints 
are not all independent, since if we add all the constraint 
equations together, we get 



=E"-'E' 



(mod m) (A8) 



which is automatically satisfied as long as 
(mod m). (This is a special case of the identity 
We conclude that the number of constraints in (IA7I) is 



1 



(A9) 



On the other hand, the number of free parameters is just 
the total number of links 



= Alink 



(AlO) 



Combining these two calculations, we see that the total 
number of configurations of n^g' satisfying (jA7[) is given 
by 



config 



,Wii„k-JV,it„ + l 



(All) 



We now substitute this into our formula (jA6|) for the 
degeneracy, obtaining 



D = 

We then note that 

Mink - 



config 



,Afli„k-Afsitc-Afplaq-|-l 



N. 



site 



iVplaq = -1 



(A12) 



(A13) 



either by direct computation or by Euler's general for- 
mula V — E + F ~ X- We conclude that D ~ 1: there is 
a unique eigenstate \qs,bp) for each configuration satis- 
fying J2s1s = (mod m). 



2. Periodic (torus) geometry 

The calculation in the torus geometry proceeds sim- 
ilarly to the open boundary condition case discussed 
above. The expression (jASp for the degeneracy is still 
applicable in this case, and we can still deduce that the 
only nonvanishing terms in this sum are the ones where 
Ajss' = (mod m) everywhere. However, unlike the 
open boundary condition case, there are now m terms 
satisfying Ajss' = (mod m) — namely the terms where 
jp is constant. Therefore, in the torus geometry case, 
[5l) reduces to 



D 







J P \ s 



r-Tr{l\sP,.) iiUpbp = l 
otherwise 



(A14) 



Just as in the open boundary condition case, we can com- 
pute Tr{Y[s Pq,) by counting the number of configura- 
tions Uss' satisfying (jA7|) . and again this number is given 
by (PTTj) . Substituting (KTH into (|M4)) . and specializ- 
ing to the case Yip bp = 1, wc find 



D = 



-N. 



config — 



A'link-iVsitc-A^pIaq + 2 



(A15) 



At this stage, there is again a difference from the open 
boundary condition case. Instead of (|A13|) . we have the 
modified relation 



iVplaq = 



(A16) 



(again, this can be derived directly or via Euler's formula 
V — E + F = x)- We conclude that there are D = de- 
generate eigenstates \qs,bp) for every configuration satis- 
fying Yip bp = 1 and qs = (mod m). In particular, 
the ground state \qs = 0,6p = 1) is m^-fold degenerate. 



Appendix B: Eigenstate degeneracy of 3D lattice 
boson model 



1. Open boundary condition geometry 

We use the equation (jASp in an almost unchanged 
form, the only difference being that the quantity Ajss' is 
now defined as jp^ — jpj + jpa —jPi , where Pi are the four 
plaquettes meeting at the edge (ss') and the labeling is 
clockwise as we go around the link (i.e. there is a relative 
minus sign between plaquettes which are perpendicular 
to each other. This is a different convention from the 2D 
case.). The only nonvanishing matrix elements are those 
for which: 



Ajss' = (mod m) 



(Bl) 



The question is thus reduced to the problem of counting 
possible assignments {jp} that satisfy those conditions. 



To count the allowed configurations of {jp}, wc will con- 
struct an explicit mapping from choices of {jp} satisfying 
Eq. ()Bip to configurations of Z^. spins living at the cen- 
ters of the cubes on the lattice (i.e. spins on the dual 
lattice). One can picture jp = as the absence of a do- 
main wall between the cubes; more generally the value 
of jp assigns a label to the face separating any pair of 
cubes, which wc will later identify with the change in the 
Zjn spin between these cubes. The conditions (jBl[) en- 
sure that domain walls can never end on an edge - and 
further, that the net change in the spin along any closed 
curve is a multiple of m, so that an allowed configuration 
of {jp} does indeed specify a set of domain walls between 
Zm spins on adjacent cubes. 

The precise mapping is as follows. Let us focus on 
some link {ss') incident at the corner of the sample ~ 
there are only two plaqucttcs P, P' bordering that link 
and Ajss' = jp — jp' , so in order for it to contribute 
to the trace we need jp = jpi which can assume m val- 
ues. We place the corresponding value of jc = jp in the 
corner cube of the lattice. We assign a value to all other 
cubes sequentially using the following algorithm: start in 
the corner cube and construct a closed directed loop go- 
ing through a number of other cubes. If C, C" are cubes 
in the loop separated by a plaquette P and the value jc 
is already assigned, then jc = jc i jp (mod m). We 
adopt the convention that we always pick up a positive 
sign as we leave the corner cube and later it alternates 
as shown in figure The consistency is guaranteed 

by the fact that as wc go around a single link {ss') and 
return to the cube C we pick up a contribution of pre- 
cisely zbAjs^/ = (mod m). More generally, if we pick 
any closed loop starting at C we can decompose it into a 
number of small loops around individual links enclosed by 
the big loop, all of which contribute (mod m), which 
can also be seen in figure (|9]). Thus the assignment of 
spins is consistent. Once the corner cube has been as- 
signed a value the mapping is unique. Conversely, given 
a configuration of {jc} in the cubes we can reconstruct 
the corresponding assignment of {jp}- Therefore there 
is a bijection between the plaquette configurations and 
the cube configurations - of which there are m^"^^" . We 
conclude that there are m^'^"'"' terms contributing to the 
sum in equation (|A5|) . 

For each plaquette configuration {jp}, we must now 
evaluate the product 

The evaluation of the trace is identical to the 2D case, 
and again gives TT^^'infe-^ff^+i. The factor Op^p' 
always equal to 1 - as long as Opgc bp = 1 for each cube 
C. (Using this relation, we can simultaneously reduce the 
labels on all of the faces of a given cube by any integer j 
(mod m). This preserves the condition (jBl[) everywhere. 
Since the allowed configurations constitute closed domain 
walls, a series of these reductions can be used to reduce all 
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FIG. 9. An example of mapping to spin configuration for m = 
3. The numbers in red correspond to a choice of {jp}- The 
base of the loop is the lower left corner, which for concreteness 
has been assigned a value 0. The sign with which we pick up 
the appropriate jp (or equivalently with which the plaquette 
contributes to Ajg^i for a specific link (ss') enclosed by the 
loop) has been coded in blue for -I- and green for — . The 
involved hnks have been colored red. 

of the labels to 0, which proves that Y[p bPp = Yip ^p = 1 
for the allowed configurations {jp})- Thus: 

The exponent can be evaluated by explicit counting or 
using the generalized Euler formula: 

A^cubc - iVplaq + Mink - A^sitc + 1 = 0. (B4) 

Thus wc have shown that D = 1, provided the constraint 
JIpj^p &p = 1 is satisfied for all cubes. 

2. Periodic (torus) geometry 

Next, we repeat the above counting argument for a 
system with periodic boundary conditions in the x, y, 
and z directions. It is instructive to start with the open 
boundary case again. We may divide the configura- 
tions {jp} satisfying (jBip into bulk and boundary parts: 

boundary 

and {jp}pebuik- Assume wc have fixed 
some bulk configuration. In the open geometry, fixing 
one of the boundary plaquettes (say, in the corner) to 
one of the m possible values automatically fixes all other 
boundary plaquettes. This is because the links {ss') on 
the boundary always have only two incident boundary 
plaquettes (and possibly a bulk one, but it is already 
fixed). Hence, the total number of configurations is m 
times the number of bulk configurations. 

Let us now define the periodic geometry by identifying 
opposite faces of the boundary. This allows us to inherit 
the notion of bulk and boundary configurations from the 
open case. Note that the two cases do not differ at all 
in the bulk - the allowed configurations satisfying (jBl[) 
are the same. However, at the boundary there are more 
possible configurations if we impose periodic boundary 
conditions. Let us again fix a bulk configuration. In 
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the periodic case choosing a value for one of the bound- 
ary plaquettes only fixes the value of the plaquettes in 
the same boundary plane. The two other perpendicular 
planes still remain unfixed as each boundary link in the 
corner has four incident boundary plaquettes. Indepen- 
dently choosing a value for one plaquette in the remaining 
planes produces a factor of im?. Hence the total number 
of configurations in periodic case is times the number 
of bulk configurations or times the total number of 
configurations in the open case. Thus we have m^cubo+a 
configurations. Note also, that in the periodic case there 
is an additional algebraic constraint satisfied by opera- 
tors Bp: if we take any boundary plane i (which is a 2D 
torus embedded in the 3D one) then: 



n Bp=i. 

P^plane{i) 



The rest of the calculation is unchanged: 



(B5) 



D = „^^cubo + 2-JVpl,q + JVli, 



-1 = m\ 



(B6) 



since in the periodic case we have A'cubc — -^piaq + Mink ~ 
-Mite = cither by explicit counting or using the Eulcr 
formula. Thus there is a D = ■mP degeneracy, provided 



the constraints ([M)) and (|B5|) are satisfied. In particular, 
the ground state \qs = 0,bp = 1) has an degeneracy. 



Substituting these expressions into the edge action (jl04p , 
we find 



L = -^d^^iiKudt^j - Vijd^^j) 
47r 

ZTT 



(C2) 



where 



/o 1 \ 




( \ 


10 




lejm 


10 


{2k + m)e/m 


yo -1/ 




\(2k + m)e/mj 



(C3) 



and V = W'^VW. In these variables, the perturbation 
becomes 

U{Ao) = 2U{x) cos($2 - a{x)) (C4) 
We next assume that that the interactions on the edge 
are tuned so that V = vSjj (we can make this assumption 
without any loss of generality since the velocity matrix 
is non-universal and can be modified by appropriate per- 
turbations at the edge). In this case, the edge theory 
can be written as a sum of two decoupled actions, one 
involving $i, $2 and one involving $3, $4: 



^34 



(C5) 



where 

L12 



(C6) 



and 



Appendix C: Relationship between the bosonic edge 
theory (|104p and the free fermion edge theory (|59|) 



In this section, we show that the perturbation U{Aq) 
(|109p where Ao = (1, 0, 0, 0) can gap out two of the edge 
modes of (|104p . leaving behind exactly the free fermion 
edge ([55]) . We accomplish this via a change of variables. 



/l \ 

1/m k/m k/m 

10 

yo 1 y 



(Cl) 



L34 = l-9.^$3(<9t$3 - vd:^^3) 

47r 

(2k + m,)e „,. „ ,~ ; . , 
2m7r 



(C7) 



It is now easy to analyze the effect of the perturbation 
U{Ao) (jC4p : this term gaps out the non-chiral Luttingcr 
liquid described by L12 by freezing the value of $2 (at 
least if U{x) is large). The resulting edge then has only 
two gapless modes, and is described by L34. On the 
other hand, it is easy to check that L34 is nothing but 
the bosonized description of the free fermion edge theory 
([59)) . with e**^,e~'*'' corresponding to the fermion cre- 
ation operators dj, o?| respectively.— >2i We conclude that 
the perturbation U{Aq) docs indeed gap out two of the 
modes of (|104p leaving the free fermion edge (|59|) . 
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